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I.  INTRODUCTION 


A.  BACKGROUND 


Microwave  circuit  elements,  such  as  semiconductor  devices,  cannot  be  connected 
directly  to  the  terminals  of  a  microwave  network  analyzer  for  measurement  of  scattering 
or  impedance  coefficients  but  must  be  coupled  through  a  chain  of  hardware  elements 
comprising  connector-adapters,  microstrip  sections,  and  bondwire  connections.  The 
basic  problem  under  consideration  is  the  removal  of  the  effects  of  intervening  hardware 
from  measurements  made  at  the  external  terminals  of  an  embedded  device  or  circuit 
component  to  obtain  the  intrinsic  device  characteristics.  Accurate  measurement  data 
on  a  device  is  essential  for  successful  circuit  design. 

Earlier  work  on  device  de-embedment  has  been  coupled  with  undesirable  restrictions 
such  as  the  assumption  of  physical  symmetry  of  hardware,  inadequate  provision  for  re¬ 
alistic  measurement  standards,  or  neglect  of  statistical  treatment  of  data  for  improve¬ 
ment  of  accuracy.  Consequently,  a  need  exists  for  a  practical  treatment  of  this 
measurement  problem  which  can  be  implemented  in  typical  laboratory  circumstances, 
with  acceptable  accuracy  of  results. 

In  order  to  measure  the  scattering  parameters  of  a  device,  the  device  is  often 
mounted  in  a  fixture.  This  fixture  has  a  set  of  scattering  parameters  associated  with  it 
so  that  the  parameters  measured  at  the  terminals  of  the  fixture  are  the  scattering  pa¬ 
rameters  of  the  device  transformed  by  the  embedding  network  of  the  fixture. 

Thus,  to  obtain  the  scattering  parameters  of  the  device  itself,  the  device  must  be 
"de-embedded"  from  the  fixture.  Techniques  for  using  a  network  analyzer  (N.A.)  to  de- 
embed  the  device  from  the  test  fixture  have  been  discussed  earlier  [Refs.  1,2,10,11]. 

As  long  as  the  device  interface  connection  is  coaxial  line  or  waveguide,  de¬ 
embedding  is  relatively  straightforward  because  standards  are  available.  But  if  the  device 
is  constructed  without  a  coaxial  or  waveguide  interface,  as  in  the  case  of  microstrip  or 
coplanar  line  systems,  de-embedding  becomes  more  difficult. 

Extracting  the  device  parameters  from  the  overall  measured  parameters  is  strongly 
dependent  on  the  ability  to  characterize  the  fixturing.  Figure  1  is  the  block  diagram  and 
signal  flow  graph  of  a  typical  N.A.-  fixturing  setup.  The  problem  is  to  find  the  best 
characterization  networks  A  and  B  so  as  to  remove  their  contributions,  yielding  the  de¬ 
vice  parameters,  S. 
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Figure  1.  Block  Diagram  and  Signal  Flow  Graph  of  a  Typical  Fixturing  Set-up 
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The  scope  of  the  proposed  research  will  be  to  develop  a  working  procedure  for  ear¬ 
ning  out  the  removal  of  embedded  device  scattering  matrices  from  network  analyzer 
data.  The  procedure  developed  will  be  supported  by  experimental  measurements  on 
typical  microwave  circuits  fabricated  on  planar  substrates. 


II.  CHARACTERIZATION  OF  THE  FIXTURES 


A.  CHARACTERIZATION  CHOICES 

Several  papers  present  characterization  techniques  that  are  in  use  at  various  labo¬ 
ratories.  There  are  four  main  choices  of  N.A.-fixturing  characterization  as  shown  Fig. 
2  (Ref.  1], 


Figure  2.  The  Four  Main  Choices  of  N.A.-Fixturing  Characterization. 


Method  (1)  combines  a  variety  of  measurements,  including  the  transmission  between 
fixtures  connected  back-to-back  and  reflection  measurements  on  the  fixtures  terminated 
in  microstrip  stubs.  In  this  way,  data  is  assembled  to  enable  the  scattering  coefficients 
of  the  embedded  device  to  be  converted  to  the  device  scattering  matrix.  Method  (1) 
therefore  contains  features  of  both  methods  (3)  and  (4)  below,  and  represents  also  an 
early  example  of  the  use  of  microstrip  stubs  as  known  reference  standards. 

In  method  (2),  the  measured  S-parameters  are  fit  to  an  assumed  circuit  model  to 
calculate  the  device  parameters. 


Method  (3)  has  the  advantage  that  only  three  simple  test  pieces  in  the  device  me¬ 
dium  are  required.  More  details  are  given  in  Chapter  V,  part  B. 

From  the  above  choices,  method  (4)  was  selected  for  investigation.  In  this  method, 
three  or  more  known  reflection  coefficients  are  measured  at  each  port.  The  reference 
load  reflection  coefficients  are  derived  from  independently  measured  or  calculated  non¬ 
zero  values.  Reference  2  used  a  chip  capacitor,  wire  bonded  to  microstrip  on  alumina, 
with  capacitance  measured  at  low  frequency  (LF).  The  overall  reflection  coefficient  was 
then  measured  versus  frequency.  To  obtain  multiple  calibration  measurements,  a  piece 
was  scribed  off  the  chip,  and  the  measurement  cycle  repeated  twfice  more  at  each  port. 

In  other  cases,  the  reference  load  was  a  high-Q  varactor  chip,  again  measured  as  a 
function  of  three  values  of  reverse  bias  at  LF.  The  chip  was  then  wire-bonded  to  the 
microstrip  conductor  and  the  overall  reflection  coefficients  of  the  fixture-chip  combina¬ 
tion  versus  frequency  were  measured  at  each  value  of  reverse  bias. 

In  the  present  work,  the  reference  load  is  a  section  of  the  microstrip  line  itself.  Three 
different  lengths  of  microstrip  line  yield  the  three  different  reflection  coefficients  for  each 
port  by  separate  measurements. 

B.  THE  TEST  FIXTURE 

The  S-parameters  of  the  test  fixture  can  be  determined  by  measuring  a  set  of  known 
impedances  in  place  of  the  embedded  device.  This  process  is  known  as  "unterminating." 
Once  the  fixture's  S-parameters  have  been  determined,  the  device  is  measured  in  the 
calibrated  test  fixture,  and  the  device  parameters  are  "de-embedded"  from  the  data  to 
give  an  accurate  picture  of  device  performance. 

The  signal  flow  diagram  and  block  diagram  in  Fig.  3  define  the  input  and  output 
characteristics  of  the  test  fixture. 

Three  sets  of  S-parameters  must  be  determined  when  an  unknown  device  is  charac¬ 
terized  as  a  two-port  network  mounted  in  a  test  fixture.  First,  the  S-parameters  of  the 
input  (A)  and  the  output  (B)  sides  of  the  fixture  must  be  evaluated  before  the  S- 
parameters  of  the  device  itself  can  be  deduced.  If  the  test  fixtures  are  known  to  be  sym¬ 
metrical,  then  only  one  measurement  is  needed  to  determine  the  S-parameters  of  the  two 
sides  of  the  fixture. 

In  the  present  case,  however,  unsymmetrical  fixtures  are  chosen,  requiring  six  dif¬ 
ferent  measurements  and  calculations  for  both  fixtures  independently.  To  solve  this 
problem,  the  input  half  of  the  fixture  is  first  terminated  at  the  device  plane  with  an  ar¬ 
bitrary  load  of  reflection  coefficient,  T^,  followed  by  other  loads  with  different  reflection 


coefficients.  Then  the  output  half  of  the  fixture  must  also  be  considered  with  same 


method. 


Figure  4.  One  Port  De-  Embedment 


For  the  one  port  case,  as  shown  above,  the  reflection  coefficient  is  the  data  needed 
since  the  impedance  can  be  calculated  from  it: 
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From  the  definition  of  the  scattering  matrix,  the  input  reflection  coefficient  is: 


b\ 
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where,  for  load  fL, 
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Thus, 


b,  =  S 


21 


1  -  $22^ 


Using  the  above,  we  get: 

■^12^21^2. 


r  n  =  s„  + 


(i  —  s22r  j 


(2.4a) 


(2.5) 


Here,  5„  can  be  found  by  measuring  a  matched  load  in  the  fixture  (TL  =  0).  If  we 
assume  that  this  measurement  gives  r,„  =  rML,  then  we  have  S„  =  TML.  Using  the  result 
in  Eq.  (2.5), 

(r,„  -  rW£j(i  -  522rj  =  s12s2.r,  (2.6) 


From  this, 


^i2^2i^l  +  s22r  L(r^  —  r  M[)  —  (r^  —  r  Ml) 


(2.7) 


Two  more  measurements  with  known  loads  will  provide  enough  data  to  determine 
Sn  and  (SUS21).  Suppose  we  measure  with  an  open-circuited  load  (T^+l).  and  a 
short-circuited  load  (Tt  »  -1).  Then, 


*>2*21  +  -^(r OC  ~  1~ML)  ~^OC~^ML 

(for 

rL=+i) 

(2-8) 
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(2.9) 

These  can  be  solved  for 
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(2.11) 


With  these  two  quantities  known,  the  equation  for  de-embedding  a  general  unknown 
can  be  written  from  Eq.  (2.7) 


rL 
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(2.12) 
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But  in  the  actual  case,  it  is  very  difficult  to  construct  the  matched  load  and  ideal 
shorted  and  open  terminations  required  by  this  procedure.  Therefore,  we  rearrange  Eq. 
(2.12)  to  matrix  form  for  three  different  known  loads  and  their  reflection  coefficients. 
From  Eq.  (2.5), 


_  ^1  i  ^12^21 

rL  ~  i-s22rL 
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Rearranging  this  gives: 
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where, 

A  ==  S11S22  ~  ^12^21- 


(2.14) 


(2-15) 


This  equation  then  involves  the  three  unknowns  S„,  Sn,  and  A  .  These  three  unknowns 
can  be  determined  by  making  three  independent  measurements  of  r,n  with  three  different 
known  loads.  In  matrix  form,  the  problem  can  be  expressed  as: 
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To  find  Su,  Sa  and  A,  solve: 
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(2.17) 


Tin  Tu  and  rL1  are  three  different  calculated  reflection  coefficients  at  the  terminals 
of  the  three  different  loads.  Also,  T,*,,  and  rM  are  three  reflection  coefficients 
measured  from  the  N.A.  with  respect  to  the  three  different  loads.  Now  the  problem  re¬ 
quires  accurate  values  of  rt„  ru,  rL3  and  r/ft„  T^.  We  consider  in  the  next  chapter 
the  computation  of  fu,  and  fi3  for  reference  microstrip  line  loads. 


i  i>i  <i  .u  ji  ji  . 


III.  MICROSTRIP  LINE  FOR  THE  FIXTURE 


A.  MICROSTRIP  LINE 

The  definition  of  an  effective  relative  dielectric  constant,  allows  us  to  treat 
microstrip  line  like  a  TEM  system.  Relations  between  line  and  freespace  parameters 
then  hold  (which  are  typical  of  TEM  systems).  With  er,  replacing  the  conventional  rel¬ 
ative  dielectric  constant, 
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where  Xm  and  vm,  are  the  wavelength  and  velocity  of  the  waves  on  the  microstrip  line  and 
and  c  are  the  free  space  values  at  the  frequency  of  operation. 

In  Eqs.  (3.1),  Lmt  and  Cmi  are  the  inductance  and  capacitance  per  unit  length  of  the 
microstrip  line.  Z„  is  its  characteristic  impedance,  and  Z^,r  is  the  value  Z0  would  have  if 
the  dielectric  layer  were  replaced  by  air.  To  the  lowest  order  of  approximation,  the  value 
of  would  be  the  average  of  the  relative  dielectric  constant  of  the  substrate  and  that 
of  air,  namely: 


e 


re  ~~ 


i  +  tT 
2 


(3.2) 


More  exact  values  have  been  obtained,  largely  through  solutions  of  the  quasistatic 
model,  corrected  to  match  measured  data. 


Figure  5.  Microstrip  Line  Dimensions. 


In  case  of  lines  with  line  width  ratio  u  a  -7-  =  0. 1  to  10, 
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If  microstrip  were  a  true  TEM  system,  the  velocity  of  waves  on  the  line,  the 
impedance  Z„  and  other  line  parameters  would  be  independent  of  frequency.  The 
quasistatic  or  zero  frequency  values  of  the  parameters  are  good  approximations  at  fre¬ 
quencies  up  to  the  order  of  a  few  gigaherz.  When  the  wavelength  becomes  comparable 
with  the  cross  sectional  dimensions  of  the  line,  vm,  and  the  line  parameters  are  no  longer 
frequency  independent.  Therefore,  dispersion  of  effective  relative  dielectric  constant  is 
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where  h  is  the  substrate  height,  f  is  frequency  and  m  is  the  permeability  of  free  space. 

B.  MICROSTRIP  STUB-RESONANCE  EVALUATION  OF  DIELECTRIC 
CONSTANT 

It  was  necessary  to  determine  the  substrate  dielectric  constant  tr  of  the  microstrip 
line  substrates  used  for  measurements  in  the  present  work. 


The  steps  are: 

1.  Measure  the  resonant  frequency  of  a  quarter  wave  length  stub.  (Suhstrate  height 
is  h=  1.5  mm,  line  width  ratio  u  =  3.0  (where  u  =  w  /  h)), 


2.  Measure  the  value  of  S2U  transmission  coefficient,  with  respect  to  frequency  using 
the  N.A., 

3.  Find  the  frequency  (f0)  at  which  maximum  dB  transmission  loss  due  to  resonance 
occurs,  as  shown  in  Fig.  7, 


Figure  7.  Determination  of  the  Max.  dB  Loss  Point. 


4.  Calculate  e„  from  Eq.  (3.1a).  It  is  normally  required  to  account  for  the  fringing 
capacitance  at  an  open  microstrip  termination.  Computation  by  standard  methods 
[Ref.  3:  p.  138]  showed  that  the  open-end  length  correction  was  less  than  1.0%  for 
stub  lengths  greater  than  50  mm.  This  had  negligible  elfect  on  the  calculated  e„ 
and  was  neglected. 

^•0  =  £re  (3-4) 
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From  the  equations  above, 

Through  the  several  measurements  as  shown  in  Table  1,  the  average  e„  ob¬ 
tained  by  calculation  has  the  value  of  2.833., 


Table  1.  THE  MEASUREMENT  OF  EFFECTIVE  RELATIVE 
DIELECTRIC  CONSTANT 


Length  (mm) 

MG  Hz) 

98.0 

0.452 

2.867 

80.0 

0.564 

2.763 

60.0 

0.746 

2.808 

40.0 

1.102 

2.895 

5.  From  stub  resonance  data,  the  corresponding  is  calculated.  Eq.  (3.3)  was  in¬ 
verted  by  trial  and  error  to  find  t,  *  3.516  from  the  experimental 

C.  CHARACTERISTIC  IMPEDANCE  OF  MICROSTRIP  LINE 

Expressions  for  the  characteristic  impedance  of  microstrip  lines  have  been  developed 
through  numerous  successive  approximations,  based  principally  on  TEM  models  having 
the  specific  microstrip  line  cross  section,  with  correction  terms  added  to  obtain  agree¬ 
ment  with  experimental  observation. 


(3-8) 


„  60  r  (0.99u+  1.795)(0.994u  +  5.80)  -| 

°"v^7  L  U(w-r  1.268) 

where  the  synthesis  formula  for  linevvidth  ratio  u  =  w,  h  is: 


»  -  f  - -J  V .<(0.W7  +  ^4^- )  +  1.32 1(  1  +  U. ) 


with 
A  = 
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exp(Z0 - gg - )-l 


(3.9) 


(3.10) 


From  Eqs.  3.3,  3.6  and  3.8,  one  can  easily  get  the  characteristic  impedance  for  each 
microstrip  section.  For  example,  for  the  given  characteristic  impedance  of  50  ft  ,  the 
microstrip  line  width  ratio  is  u  =  2.283  and  the  width  is  3.43  mm.  In  the  construction 
of  microstrip  stub  loads,  it  is  necessary  to  choose  high  characteristic  impedance  for  ac¬ 
curate  measurement,  since  high  impedance  leads  to  narrow  linewidth  and  a  well-defined 
point  of  attachment  of  the  stub.  Construction  problems  arise  in  making  microstrip  line 
with  very  small  width.  For  example,  if  we  want  to  make  a  microstrip  line  stub  of  90  ft 
characteristic  impedance,  the  required  width  is  about  1.2  mm  for  given  h  =  1.5  mm.  It 
is  difficult  to  handle  the  adhesive  foil  material  used  for  the  microstrip  itself.  For  this 
reason,  a  compromise  value  was  chosen,  using  2.0  mm  width  of  microstrip  stub  for  given 
substrate  height  h  =  1.5  mm.  Then, 

h  =  1.5mm 

w  —  2.0mm 

also  u  =  1.33. 

h 


Therefore, 


Z0s  =  68.2ft, 


which  was  chosen  as  a  working  value  for  stub  impedance. 
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MICROSTRIP  UNE 

CHARACTERISTIC  IMPEDANCE  :  50  Ohm 

WIDTH  :  3.43  mm 

OPEN  STUB 

CHARACTERISTIC  IMPEDANCE  :  68.2  Ohm 

WIDTH :  2.00  mm 

connector 


50  ohm 


68.2  ohm 

2.0  mm 


Figure  8.  Experimental  Setup  for  Characterization  of  the  Fixture 


D.  STUB  LOAD  IMPEDANCE  AND  REFLECTION  COEFFICIENT 

Consider  the  transmission- line  impedance- transformation  law: 


-In  -t-Ol 


["  ZL  cos  pi  +jZ0S  sin  pi  1 
|_  Zqs  cos  fil  +  ZL  sin  pi  j 


(3-11) 


From  Eq.  (3.11),  when  ZL  goes  to  oo  (open  ended),  the  input  impedance  of  an  open  - 
ended  stub  becomes: 
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(3.15) 


Once  the  Z0l  is  chosen,  we  can  determine  the  load  impedance  of  three  different  loads. 
For  example,  given  Z„,  =  68.2  ohms  and  microstrip  line  lengths  of  4  cm,  3  cm,  2  cm,  the 
load  impedances  are  tabulated  as  shown  Table  2. 


Table  2.  THE  LOAD  IMPEDANCE  OF  REFERENCE  STUB  LOADS 


Frequency(GHz) 

L  ■■  4  cm  ( ZLl ) 

1.00 

-j  11.06 

-j  38.44 

-j  80.15 

2.00 

j  204.8 

j  41.28 

-j  11.06 

3.00 

-j  35.69 

-j  1230 

j  41.28 

4.00 

j  91.05 

-j  35.69 

j  204.81 

5.00 

-j  70.73 

j  44.23 

-j  168.99 
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The  three  different  reflection  coefficients  of  the  given  arbitarv  loads  are: 


ru- 


Zjj  ~ 

Zu  +  Z0 


(3.16) 


where  i  =  1,  2,  3.  Therefore,  from  the  given  stub  impedances  and  Z0  =  50Q,  the  re¬ 
flection  coefficients  are  tabulated  in  Table  3. 


Table  3.  THE  STUB  REFLECTION  COEFFICIENTS 


Freq.(GHz) 

1.00 

-0.91  -  j0.42 

-0.26  -  j0.97 

0.44  -  j0.90 

2.00 

0.S9  +  j0.46 

-0.19  + jO.98 

-0.91  -  j0.42 

3.00 

-0.32  -  j0.95 

1.00  -  j0.05 

-0.19  +  j0.98 

4.00 

0.54  +  j0.84 

-0.32  -  jO.95 

0.89  +  j0.46 

5.00 

0.33  -  j0.94 

-0.12  + j0.99 

0.84  -  j0.54 
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IV.  APPLICATION  OF  THE  DE-EMBEDDING  PROCEDURE 

A.  MEASUREMENT  OF  THE  EMBEDDING  FIXTURE 

The  reflection  coefficients  required  are  measured  by  the  HP  8409B  vector  network- 
analyzer.  One  can  find  the  reflection  coefficient,  r,n,  of  the  stub  terminated  fixtures  by 
measurement  of  the  scattering  matrix  element  S„  using  the  network-analyzer.  The  out¬ 
puts  of  the  N.A.  represent  the  magnitude  (dB  loss)  and  phase  angle  (degrees)  of  S„. 
Therefore,  the  outputs  of  the  N.A.  must  be  changed  to  complex  notation  using  the 
simple  modification: 

ri>t  =  (10)  20  cos(argSM)+y(10)  20  sin(argSu).  (4.1) 

The  outputs  of  the  N.A.  and  are  obtained  for  each  load  on  one  fixture  (A)  as  shown 
in  Figs.  9-a  to  9-c.  When  we  consider  the  other  fixture  (B),  all  steps  are  exactly  same 
as  for  fixture  A  except  the  fixture  shape  itself,  shown  in  Figs.  10-a  to  10-c.  Tables  for 
fixture  A  and  fixture  B  data  are  available  in  Appendix  A. 

B.  S-PARAMETERS  OF  THE  FIXTURE 

To  consider  the  de-embedding  calculation,  we  return  to  the  original  matrix  equation. 
The  success  of  this  characterization  technique  depends  heavily  on  the  measurement  of 
the  load  impedances  used  in  the  unterminating  process  for  the  fixtures.  Once  an  input 
reflection  coefficient  has  been  measured  with  three  different  impedance  standards,  the 
S-parameters  of  the  test  fixture  can  be  computed  using  the  matrix  introduced  earlier. 
A  word  of  caution  must  be  offered  about  the  convergence  of  this  apparently  simple 
problem. 

The  matrix  that  must  be  inverted  consists  of  elements  very  near  unity,  a  situation 
often  used  to  illustrate  some  of  the  numerical  inversion  methods.  In  this  case,  the 
Gauss-Jordan  method  with  full  pivoting  is  recommended.  Pivoting  improves  the  stabil¬ 
ity  of  the  numerical  solution  by  interchanging  rows  and  columns  in  the  Gaussian  elim¬ 
ination.  In  the  notation  of  Eq.  (2.17),  the  scattering  parameter  S0  is  changed  to  A,,  and 
Bij  for  fixture  A  and  B  (respectively): 
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Figure  9.  Measurement  of  the  Load  for  Fixture  A 
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Figure  10.  Measurement  of  the  Load  for  Fixture  B 
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where 

A  =  AUA22  ~  ^12^21* 


rZ.l 
r  ini 
^L2 
Hw3 


(4-2) 


(4.3) 


The  computed  values  for  Alx  and  An,  combined  with  a  knowledge  of  A  ,  can  be  used  to 
determine  the  AX2  A2l  product.  Note  that  the  fixture  is  reciprocal,  so  An  =  A2X  and  for  the 
other  fixture,  Bl2  =  Blv  Thus,  the  magnitudes  of  these  terms  can  be  calculated  by  taking 
the  square  root  of  the  product  term  in  the  expression  for  A  . 

The  question  of  phase  is  left  open,  but  is  not  essential.  Since  the  Ai2,  A 2I,  Bl}  and 
B2l  terms  always  appear  in  product  form,  it  does  not  matter  which  of  the  two  possible 
phase  angles  one  chooses,  as  long  as  the  phase  angles  for  A2I  and  B2l  are  approximately 
the  same.  All  of  these  calculations  are  shown  in  Appendix  E  using  the  Fortran  program 
shown.  This  program  combines  the  measurement  and  calculated  values,  that  is,  the 
outputs  of  the  N.A.  are  used  in  the  calculation  as  a  form  of  datafile.  The  resulting  S- 
parameters  of  these  fixtures  are  shown  in  Appendix  B. 


C.  DE-EMBEDDING  PROCEDURE 
1.  Background 

When  a  two-port  network  is  connected  between  fixture  A  and  B,  the  signal  flow 
graph  for  the  system  may  be  drawn  as  shown  in  Fig.  11.  Parameters  AB  and  BA  are 
mainly  contributed  by  the  finite  isolation  of  the  internal  switches  of  the  network  ana¬ 
lyzer.  The  terms  S,„  S3„  Sia,  and  S22  are  the  scattering  parameters  of  the  device  under 
test.  The  fixture  parameters,  Au,  A2l,  A, 2  and  Au,  can  be  determined  by  reflection 
measurements  at  port  I,  as  discussed  earlier.  Parameters  Ba,  Bu  and  BnB2i  can  be  eval¬ 
uated  by  similar  measurements  at  port  2.  Two  port  transmission  measurements  with  the 
device  under  test  replaced  by  a  through-line  are  used  for  evaluating  parameters  AB  and 
BA,  and  products  A2lB2l  and  AnBn.  When  a  through-line  of  length  /,  having  zero  re¬ 
flections  at  both  ends  and  a  transmission  T  =  e"*',  is  connected  between  the  two  ports, 
the  signal-flow  graph  of  Fig.  11  can  be  simplified  to  the  one  shown  in  Fig.  12. 
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Figure  11.  Signal  Flow  Graph  for  Measurements  on  a  Two  -  Port  Network 


The  measured  transmission  Tm  can  be  expressed  as 


Tm-  BA  + 


I  -AaBuT‘ 


(4.4) 


When  the  error  parameters  Aa  and  Bn  have  been  determined  by  the  reflection  measure¬ 
ments,  two  transmission  measurements  must  be  performed  for  the  determination  of  BA 
and  The  value  of  BA  can  be  found  by  making  a  transmission  measurement  when 
no  connection  exists  between  the  ports  (that  is  T  =  0).  The  parameter  AlxBn  can  be 
determined  from  the  measurement  of  a  through  -  line  using 

Ai\Bi\  = - — -  (4.5) 


where  T  =  e-w,  \{i  being  the  propagation  constant  of  the  through-line.  The  error  pa¬ 
rameters  AB  and  Al2B,3  can  be  determined  by  making  a  similar  measurement. 

Now  that  the  scattering  characteristics  of  the  embedding  networks  have  been 
measured,  the  scattering  parameters  Mn,  Mu,  Mn  and  A/2,  can  be  corrected  for  the  em- 
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Figure  12.  Signal  Flow  Graph  with  a  Through  -  Line 

bedding  networks  to  obtain  S,„  S2,  Su  and  S2l.  Elements  Mtj  represent  the  scattering 
parameters  of  both  fixtures  and  device,  namely: 


^11^12^21(1  ”  ^22^11)  +  (^21^12^11^12^21) 
Mn  =  An  + - 7t - 


,,  „  i  , 

M2\  =  BA  + - jr— 


^22^21^12(1  —  ^11^22)  +  (^21^12/l22^12^2l) 
Mu  =  Bn  + - - - 

w  in.  ^12^12^12 

Ml2  =  AB  + - ^ - 

where 

D  =  1  —  Su/l22  “  ^22^11  -  *^1 2^21  ^22^1 1  +  ^22^1 1^22^1 1 


(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 
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These  equations  must  be  inverted  to  extract  the  scattering  parameters  StJ.  An  explicit 
solution  for  the  scattering  parameters  of  the  device  under  test  has  been  presented  in  Ref. 
[4] :  p.  313  .  The  solution  is  as  follows: 


G(l  +  BUH)-BUEF 

-  ,y 

(4.11) 

^  -  .V 

(4.12) 

5  __L 

521  ■  .V 

(4.13) 

H(l+AUG)-AUEF 
*22  ~  N 

(4.14) 

where 

jV=  (l  +  A22G){1  +  BUH)  -  AUBUEF 

(4.15) 

n  Wn-AB) 

(AuBn) 

(4.16) 

«  (Afa.  -  BA) 

(/42,52i) 

(4.17) 

1  ^ 

1  rM 

v=  3 
— 

n 

(4.18) 

j,  W22  “  B22) 

(BnB2x) 

(4.19) 

In  the  present  measurement,  inclusion  of  the  leakage  paths  BA  and  AB  is  not  necessary 

because  of  the  one-port  calibration  procedure  used. 

Therefore,  Eqs.  (4.16,  4.17)  are: 

£  iV/'2 

AuBn 

(4.16a) 

F. 

^21^21 

(4.17a) 

Some  variations  of  the  above  calibration  procedure  are  available  in  Ref.  [4] :  p.  315  . 


The  de-embedding  procedure  used  in  this  work  is  based  on  Eqs.  (4.11  to  4.15) 
and  Eqs.  (4.16a,  4.17a,  4.18  and  4.19)  which  use  the  9  key  parameters: 

Sn,  S|2.  S21,  S22,  E,  F,  G,  //,  iV 

and  the  foregoing  derivations.  The  complete  block  diagram  is  shown  in  Fig.  14.  If  we 
know  Av  and  B,r  the  scattering  parameters  of  fixtures  A  and  B,  and  the  parameters  of 
both  fixtures  plus  device,  MtJ,  then  the  scattering  parameters  of  the  device  itself  are  ob¬ 
tained  by  using  Eqs.  4.1 1  to  4.14.  In  order  to  test  the  de-embedding  procedure,  consider 
the  set  up  of  Fig.  13. 


Figure  13.  Device  and  Fixture  Assembly  for  the  51  Ohm  Resistor 


It  is  necessary  to  choose  a  device  that  has  known  scattering  parameters;  then  it  will  be 
possible  to  compare  the  output  of  the  de-embedding  experiment  with  the  theoretical 
parameters  of  device  itself.  The  initial  test  object  is  a  series-connected  51  Ohm  resistor 
which  has  a  small  inductance  so  that  the  impedance  of  the  51  Ohm  device  is: 


Z  —  R  +  ja)L 
▼  51  +j2n/L 


(4.20) 


where  we  may  assume  the  parasitic  inductance  associated  with  the  resistor  to  be 
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Fixture  A 


Fixture  B 


Figure  14.  The  Block  Diagram  of  the  De-Embedding  Procedure 
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The  scattering  matrix  of  an  impedance  Z  connected  in  series  with  the  trans¬ 
mission  line  is : 


(4-21) 


and  Z„  is  the  characteristic  impedance  (50  Ohm)  of  the  connected  line.  In  this  case,  the 
TOUCHSTONE  CAD  program  is  useful  in  calculating  the  scattering  parameters. 

These  theoretical  outputs  are  shown  in  Appendix  C.  Also  the  magnitude  plots 
of  the  calculated  scattering  parameters  of  this  device  are  shown  in  Fig.  15. 
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Theoretical  Outputs 
Magnitude  plots 
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Figure  15.  Theoretical  Scattering  Parameters  of  51  Ohm  Resistor 
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2.  Device  Embedment 

The  S-parameters  of  both  fixtures  plus  device  are  measured  by  the  network  an¬ 
alyzer  at  planes  AA  and  BB’  as  shown  in  Fig.  16.  Typical  measured  magnitude  data  is 
shown  in  Fig.  17  with  data  tables  in  Appendix  C. 


measurement  plane 


Figure  16.  Signal  Flow  graph  for  the  Measurement  of  the  (Fixture  +  Device)  S- 
parameters 
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Fixtures  +  51  Ohm  Resistor  Outputs 
Magnitude  plots 


Freq.  (GHz) 


Figure  17.  The  Magnitude  Plot  of  Both  Microstrip  Fixture  and  51  Ohm  Resistor 
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3.  De-embedding  the  Device 

Now  we  are  ready  for  calculation  of  the  S-parameters  of  the  device.  The  S- 
parameters  of  fixture  A,  fixture  B  and  both  fixtures  plus  device  constitute  the  necessary 
data.  Once  the  microstrip  circuit  board  is  prepared,  then  the  network  analyzer  provides 
the  three  different  input  reflection  coefficients  r,nl,  F,^  and  r;n3  with  respect  to  three  dif¬ 
ferent  calibration  loads.  Also  computations  are  carried  out  to  find  the  load  reflection 
coefficients  of  these  three  different  loads,  ru,  TL1  and  rt3.  After  following  the  same 
procedure  for  both  fixture  A  and  fixture  B,  we  can  get  the  respective  S-parameters. 
From  the  matrix  equation  (Eq.  2.17),  the  three  unknowns  Sn,  Sn  and  A  are  solved.  For 
de-embedding  the  device  itseif,  we  need  the  additional  S-parameters  for  the  device  plus 
fixture.  These  three  sets  of  S-parameters  are  utilized  to  solve  the  de-embedding  problem. 

Equations  4.11  to  4.14  are  now  ready  for  use.  Through  these  equations,  the 
de-embedded  outputs  are  obtained,  as  shown  in  Fig.  18  (the  table  of  de-embedded  S- 
parameters  is  shown  in  Appendix  C).  In  addition  to  the  plot  of  the  de-embedded  mag¬ 
nitudes,  the  comparison  of  the  scattering  parameters  S,„  Sn  and  Su  of  the  de-embedded 
output  with  the  parameters  of  the  embedded  resistor  and  the  theoretical  values  are  also 
available  in  Fig.  19  to  Fig.  21  (data  tables  in  Appendix  C), 
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Figure  19.  Comparison  of  the  Sll  of  the  51  Ohm  Resistor  with  the  Theoretical 
Values 
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Comparison  of  Output  S12 

(1) .  Fixtures  +  51  Ohm  Resistor 

(2) .  Theoretical  of  Resistor 

(3) .  De-embedded  output 


Figure  20.  Comparison  of  the  S12  of  the  51  Ohm  Resistor  with  the  Theoretical 
Values 
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Figure  21.  Comparison  of  the  S22  of  the  51  Ohm  Resistor  with  the  Theoretical 
Values 


V.  FURTHER  DE  -  EMBEDDING  METHODS 


A.  MEASUREMENT  OF  DEVICE  PARAMETERS  USING  A  SYMMETRIC 
FIXTURE. 

I.  Background 

As  discussed  earlier,  if  the  device  is  constructed  without  a  coaxial  or  waveguide 
interface,  de-embedding  becomes  more  difficult.  In  some  cases,  it  is  possible  to  proceed 
using  an  auxiliary  standard  whose  properties  are  calculated  or  measured  at  low  fre¬ 
quencies,  but,  in  other  cases,  no  suitable  auxiliary  standard  is  available.  In  either  case, 
if  the  fixture  is  symmetric,  a  de-embedding  procedure  that  does  not  require  the  use  of 
an  auxiliary  standard  can  be  developed  [Ref.  5]. 

The  concept  for  symmetric  de-embedding  is  illustrated  in  Fig.  22.  Fig.  22-a 
shows  a  device  inserted  in  a  symmetric  test  fixture.  We  assume  that  the  device  has 
complementary  interfaces  on  the  input  and  output  sides  so  that  when  the  device  is  re¬ 
moved  from  the  fixture,  the  two  halves  of  the  fixture  can  be  connected,  as  shown  in  Fig. 
22-b.  The  scattering  parameters  of  each  half  of  the  test  fixture  can  be  extracted  from  the 
scattering  parameters  of  the  whole  test  fixture  connected,  as  illustrated  in  Fig.  22-b. 
Then  the  device  is  inserted  into  the  test  fixture,  as  illustrated  in  Fig.  22-a,  and  the  scat¬ 
tering  parameters  are  again  measured. 

The  method  for  implementing  the  concept  outlined  above  is  as  follows.  If  we 
have  a  fixture  whose  T-parameters  in  a  50  ohm  system  are  ZTF],  then  if  \_HTF~\  re¬ 
presents  the  square  root  of  [TF], 

C  HTFllHTn  =  ZTF]-  (5.1) 

Now  suppose  we  want  to  measure  the  T-parameters  of  a  device  in  a  50  ohm  system 
CrSX  If  we  separate  the  fixture  into  two  symmetric  parts  and  inserted  the  device  to  be 
measured,  then  for  the  combination  we  would  measure  a  set  of  T-parameters  C7TS] 
given  by: 

ITFS1  =  tHTnZTSIZHm  (5.2) 


To  find  CFS[],  we  would  simply  invert  above  equation  to  get 

[75]  =  ZHTFY'ZTFS~\ZHTFYX  ■ 


(5-3) 


Figure  22.  Device  Embedded  Between  Symmetric  Fixture  and  Test  Fixture 
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This  procedure  de-embeds  the  device  parameters  [TS3  from  the  measured  parameters 
C  TFS2  and  C  TF].  The  steps  required  to  find  C  FS]  in  a  symmetric  fixture  are  as  follows: 


1.  Measure  the  scattering  parameters  of  both  fixture  plus  device,  QSFS],  using  the 
network  analyzer, 


CSFSJ - 


SFSU 

SFS2l 


SFS 

SFS 


'12 
22  _ 


(5.4) 


then  transform  to  the  T-parameters: 


[FFSIl 


-Mi 

SFS2\ 

—sfs22 

~SFS^~ 


SFS 


n 


sfs2X 

l 


SFS 


21 


(5.5) 


where  A/,  =  SFSl{SFSa  +  SFSl2SFS2l. 

2.  Measure  the  scattering  parameters  of  both  fixtures  C SF].  In  practice,  both  fixtures 
are  measured  with  the  network  analyzer  individually: 


IHSF1  =  C//SR]  =  [//5£2]  = 


IHTF2  = 


HSF- 21 
-HSFn 


HSF2X  nc>r2\ 

where  \t2  =  HSFnHSFu  +  HSFl2HSF2{ 
3.  Invert  [ HTFJ: 

HTFn 


ZHTFY'  = 


where  A/3  =  HTFnHTFu  -  HTFnHTFn. 
4.  Find  device  T-parameters  CTS]: 


crs3  = 


HSFu  HSFn 
IiSF2l  HSF22 


Then,  transforming  to  the  T-parameters: 
-Ft 2  HSFU 


HSF2l 

1 

HSF , 


A/3 

~HTFn 

A/i 


—HTF- 

A/, 


21 


//7T, 

FU 


~HTFU  HTFl2 

-1 

‘rrc„  tfsx{ 

htfxx  htfx2 

htf2X  htf22 

TFS»  TFSn 

htf2X  htf22 

-1 


5.  Transform  [TS]  to  S-parameters  CSS]: 

5,2 1 
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(5.6) 


(5.7) 


(5.8) 


(5.9) 


(5.10) 

(5.11) 

(5.12) 

(5.13) 


2.  De- Embedding  of  a  Filter  Unit 

As  a  test  demonstration  of  the  symmetrical  de-embedment  procedure,  an  8436A 
bandpass  filter  (8-12  GHz)  is  used  as  the  test  device.  All  components  of  the  network 
had  type  APC-7  test  port  connectors.  The  outputs  and  [//SP1  of  the  network 

analyzer  are  shown  in  Appendix  D.  To  get  the  de-embedded  outputs,  the  Fortran  pro¬ 
gram  shown  in  Appendix  F  is  used.  The  magnitude  plot  ol  the  filter  embedded,  the  filter 
itself,  and  the  de-embedded  filter  are  shown  in  Figs.  23  to  25.  Data  Tables  are  shown 
in  Appendix  D.  The  comparison  of  the  scattering  parameters  Su,  S12  and  S22  of  the  filter 
itself,  and  embedded  filter  are  shown  in  Figs.  26  to  28.  Tabulated  data  is  shown  in  Ap¬ 
pendix  D. 
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Figure  23.  The  S-parameter  Magnitude  Plot  of  Both  Fixture  and  Filter 
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Filter  without  Fixture  Outputs 
Magnitude  Plots 
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Figure  24.  The  S-paranieter  Magnitude  Plot  of  Filter  Without  Fixtures 
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Comparison  of  Output  S11 

(1) .  Fixtures  +  Fitter 

(2) .  Filter  without  Fixtures 

(3) .  De-embedded  output 
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Figure  26.  Comparison  of  the  Sll  of  De-embedded  Filter  with  Filter  Without  Fix- 
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Comparison  of  the  S12  of  De-embcdded  Filter  with  Filter  Without  Fix¬ 
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Comparison  of  Output  S22 

(1) .  Fixtures  +  Filter 

(2) .  Filter  without  Fixtures 

(3) .  De-embedded  output 
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Figure  28. 


Comparison  of  the  S22  of  De-enil>edded  Filter  with  Filter  N'  ithout  Fix- 


B.  TSD  (THROUGH,  SHORT,  DELAY)  METHOD 

This  procedure,  developed  by  R.  Speciale  [Refs.  6.  7  ],  is  intended  for  system  calibration 
and  systematic  error  remora  in  automated  S-parameter  measurement.  This  procedure 
uses  only  three  non*critical  reference  standards,  called  the  "Through",  the  "Short"  and 
the  "Delay",  and  leads  to  a  set  of  equations  for  which  a  fast,  non-iterative,  explicit  sol¬ 
ution  is  known  [Ref.  6]. 

1.  The  Basic  Concept 

This  procedure  uses  only  three  mechanically  simple  reference  standards,  known 
respectively  as  the  "Through",  the  "Short",  and  the  "Delay". 

The  "Through"  is  defined  as  the  direct  mutual  connection  of  fixtures  A  and  B, 
as  shown  in  Fig.  29,  when  the  uncalibrated  S-parameter  data  for  the  device  are  acquired. 
The  S-parameter  data  for  "Through"  are  therefore  acquired  with  the  measurement  lines 
to  the  device  directly  interconnected  and  the  device  removed  from  the  signal  path. 

The  "Short"  is  defined  as  a  reference  object  presenting  perfect  shorts  at  the  im¬ 
mediate  interfaces  where  the  unknown  device  is  connected  for  measurement. 

Finally,  the  "Delay"  is  defined  as  any  arbitrary  length  of  50  ohm  line  inserted 
between  the  interfaces  where  the  unknown  device  is  connected  for  measurement.  Sum¬ 
marizing,  for  the  "Through"  and  the  "Delay",  the  signal  path  is  continuous  between  ports 
1  and  2  of  the  S-parameter  unit,  through  the  couplers  used  for  connecting  the  object. 
A  direct  connection  represents  the  "Through"  and  a  finite  arbitrary  length  line  for  the 
"Delay",  in  place  of  the  unknown  object.  On  the  contrary  the  signal  path  is  blocked  for 
the  "Short"  by  perfect  shorts  appearing  at  the  interfaces  where  the  unknown  device  is 
normally  connected. 

Three  S-parameter  matrices  [SrX  and  are  acquired  at  each  desired 
frequency  during  system  calibration.  These  matrices  describe  the  "Through",  the  "Short" 
and  the  "Delay"  respectively.  Because  of  the  device  imperfections,  these  matrices  do  not 
represent  theoretically  perfect  "Through",  "Short"  and  "Delay"  elements,  which  are  as¬ 
sumed  to  exist  in  the  actual  measurement. 

Now  consider  the  error  boxes  A  and  B.  Error  box  A  is  assumed  to  be  connected 
on  the  side  of  port  1,  while  error  box  B  is  assumed  to  be  on  the  side  of  port  2,  as  shown 
in  Fig.  29. 


Figure  29.  Block  Diagram  of  the  TSD  Calibration  Procedure 


From  Fig.  29,  the  error  boxes  A  and  B  are  always  assumed  to  be  inserted  at  the 
two  ends  of  the  measurement  50  ohm  lines  connecting  the  device  to  ports  1  and  2  of  the 
S-parameter  unit.  An  algorithm  has  been  derived  wlr.h  performs  an  explicit,  non¬ 
iterative  calculation  of  the  S-parameter  matrices  [S*]  and  CSs!L  which  characterize  the 


virtual  error  boxes  A  and  B. 


2.  The  Equations  and  their  Explicit  solution 

As  discussed  earlier,  the  T-matrices  can  be  multiplied  simply  with  cascaded  ele¬ 
ments.  Therefore,  consider  Fig.  30,  where  the  T-matrix  for  the  connection  shown  is: 


(5.14) 


From  the  above,  we  get: 


L7-0]  =  cr,r'  cw3  cts]-1 


(5.15) 


Figure  30.  Device  Embedded  Setup 


where  [  TD]  is  the  corrected  T-matrix  and  Cr^os3  is  the  uncalibrated  T-matrix  For  the 
device. 

The  T-parameters  of  Fixtures  plus  device,  [r4DS],  are  given  by  the  network 
analyzer.  But  the  problem  is  to  obtain  the  T-parameters  of  each  Fixture,  [  TA 1  and 
To  do  this,  we  consider  the  error  boxes  A  and  B  as  showm  in  Fig.  31. 

From  a  formal  point  of  view,  both  the  "Through"  and  the  "Delay"  can  be  con¬ 
sidered  as  different,  Finite  lengths  of  50  ohm  lines,  cascaded  between  the  error  boxes  A 
and  B.  With  CTc,]  and  [  TD1~\,  the  T-matrix  derived  From  the  S-parameter  data  obtained 
in  these  two  calibrations,  and  and  CT*],  the  unknown  T-parameters  of  the  error 
boxes  A  and  B,  we  can  write: 


ITDJ  =  ITAHTUKTB1, 

ITD{]  =  ITA-}ITL2HTB1. 

For  the  50  ohm  lines  T,  and  L2  we  have: 


where  p,  =  r,  +  j9t. 
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(5.16) 

(5.17) 

(5.18) 


Figure  31.  Error  Boxes  of  A  and  B  with  Delay 
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where  p,  =  r2  +  j02. 

Equations  (5.16)  and  (5.17)  together  with  the  two  equations  for  the  "Short" 
characterize  boxes  A  and  B.  Also,  consider  the  situation  for  the  "Short",  as  shown  in 
Fig.  32,  we  have 
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is 


(5-21) 


where  S115  and  are  Sn  and  of  the  "Short",  and  SiiA,  S,lt  are  S„  of  box  A  and  box 
B  (det  means  matrix  determinant),  where  i=  1,  2. 

From  Eqs.  (5.16)  and  (5.17),  we  get: 
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Figure  32.  Short  Situation  of  Two  Error  Boxes 
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Now  substitute  Eqs.  (5.22)  and  (5.22a)  into  (5.20)  to  obtain: 
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(5.24) 
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and  also, 

•^12/4^21/1  =  ^11/4^22/4  ~  ^et 


(5.25) 
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Substituting  Eqs.  (5.23)  and  (5.23a)  into  (5.21), 
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where 


(5.28) 


(5.29) 


Hn 
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2  sinh  p 


(5.33) 


Here  the  quantities  Htj  and  Ktj  are  the  elements  of  the  matrices  C#]  and  CAT],  defined 


C//]  =  crD1]  itD2t 


(5.34) 


ito  =  iTD2y'  crD1] 


(5.35) 


and  p  is  defined  by: 


P  =  cosh-‘[  j-(Hn  +  //22)]  =  cosir^y(^tl  +  A:22)]  (5.36) 

Once  the  S-parameters  of  the  error  boxes  A  and  B  are  known  as  given  by  Eqs.  (5.22) 
through  (5.27),  use  the  general  conversion  formulas, 
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where  det  S  =  SUS,2  —  S,2S21. 

As  discussed  before,  the  following  calculation  to  de-embed  the  device  matrix, 

itdi  =  itay'  c tadbi  ra-1. 


can  be  carried  out  using  Eqs.  (5.22)  through  (5.27). 

3.  Discussion 

This  calibration  procedure  uses  only  three  reference  standards,  which  is  the 
minimum  necessary'  for  characterizing  the  error  boxes  A  and  B.  By  using  two  delays  of 
different  length,  for  a  total  of  4  standards,  a  check  of  intrinsic  consistency  may  be  ob¬ 
tained.  The  standards  used  are  simple  enough  to  be  inserted  at  the  location  of  the  ob¬ 
ject,  even  if  this  is  an  active  device  contained  in  some  kind  of  measurement  fixture.  In 
this  way,  all  system  errors  are  removed,  including  those  introduced  by  the  fixture. 

The  position  of  the  reference  plane  is  immaterial,  as  long  as  it  is  kept  constant 
during  calibration  and  measurement.  Also,  the  given  solution  of  the  relevant  network 
equations  is  an  exact,  explicit  solution.  Finally,  the  procedure  is  essentially  independent 
of  the  specific  system  configuration. 
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VI.  CONCLUSION 


Various  sources  have  been  cited  for  the  extraction  of  the  scattering  parameters  of  a 
device  by  removal  of  the  error  and  contributions  of  the  fixturing  setup  [Refs.  8,9,10,1 1]. 
In  Chapters  II  and  III,  the  de-embedding  procedure  using  microstrip  line  within  the 
given  conditions  was  discussed.  This  project  included  various  aspects  such  as  the 
microstrip  line  fixturing  setup  and  the  construction  and  handling  of  small  experimental 
boards.  But  in  spite  of  the  limitations  of  the  early-model  Network  analyzer  HP  8409B 
which  was  available  for  this  work,  the  best  possible  results  were  obtained  using  repeated 
measurements  on  the  same  experimental  board  to  obtain  mean  values. 

To  extract  the  scattering  parameters  of  the  device  51  ohm  resistor,  the  nine  basic 
parameters  as  discussed  earlier  in  Chapter  IV  (Eqs.  4-11  to  4-19)  were  used.  These 
equations  required  the  scattering  parameters  of  both  fixtures  since  they  are  not  sym¬ 
metrical.  That  is,  individual  measurements  and  calculations  were  required  for  both  fix¬ 
tures.  The  matrix  equation  (Eq.  2-17)  involved  the  three  unknown  scattering  parameters 
Su,  and  A  s  (S„S22  —  Sl2S21)  for  fixture  A  and  B.  That  meant  each  fixture  required 
three  equations  for  the  matrix  solution.  Three  different  measurements  and  calculations 
of  microstrip  line  wer*  used  to  set  up  the  above  matrix  entries,  i.e,  three  different  meas¬ 
ured  reflection  coefficients  from  the  network  analyzer  with  respect  to  the  three  different 
microstrip  loads,  respectively  Tml,  T(nl  and  r„,3  and  three  different  calculated  reflection 
coefficients  at  the  terminals  of  the  three  different  loads  ru,  Ti2  and  T0  were  required. 
The  computed  value  for  Su  and  S22  are  combined  with  a  knowledge  of  A  which  can  be 
used  to  determine  the  S,2S21  product.  Since  S12,  and  S21  for  the  fixtures  always  appear  in 
product  form,  the  magnitude  of  these  terms  can  be  calculated  as  the  square  root  of  the 
product  term  in  the  expression  for  A  . 

The  de-embedded  scattering  parameters  of  the  device  (51  ohm  resistor)  were  found 
through  experiment  to  be  close  to  the  theoretically  calculated  values,  shown  in  magni¬ 
tude  plots  in  Figs.  19,  20,  and  21  and  as  tabulated  in  Table  4.  As  may  be  seen  in  Figs. 
19  to  21,  the  theoretically  calculated  scattering  coefficients  of  the  51 -Ohm  resistor  device 
showed  only  a  slight  continuous  variation  over  the  entire  frequency  range  of  interest. 
The  measured  values  of  the  embedded  device  showed,  however,  rapid  irregular  fluctu¬ 
ations  with  frequency,  reflecting  the  probable  occurrence  of  internal  resonances  within 


parasitic  discontinuities  of  the  coupling  circuitry.  It  is  the  task  of  the  de-embedding 
procedure  to  remove  these  effects  from  the  measurement.  The  curves  numbered  (3)  in 
the  above  Figures  show  the  de-embedded  measurements,  in  comparison  with  the  the¬ 
oretical  curves.  As  may  be  seen,  the  de-embedded  results  followed  the  expected  values 
closely,  showing  fluctuations  which  are  assumed  to  be  due  to  measurement  error,  and 
network-analyzer  calibration  errors. 


Table  4.  THE  COMPARISON  OF  THE  OUTPUT  FOR  THE  51 
OHM  RESISTOR  IN  DB  UNITS 


Freq.(MHz) 

Theoretical 

De- 

embedded 

Embedded 

6.11 

6.35 

11.14 

2300.0 

5.54 

5.12 

10.18 

2600.0 

5.02 

4.22 

7.13 

2900.0 

4.55 

3.87 

6.21 

3200.0 

4.13 

3.19 

7.33 

3500.0 

3.76 

3.16 

11.29 

3SOO.O 

3.42 

3.91 

. 

11.59 

4100.0 

3.12 

2.66 

6.27 

4400.0 

2.86 

4.11 

6.07 

4700.0 

2.62 

2.25 

9.91 

5000.0 

2.41 

2.11 

11.95 

S,j  and  S„  terms  for  fixture  A  and  B  always  appears  in  product  form,  considering 
their  phase  is  not  important.  As  long  as  the  phase  angles  of  S2,  for  fixture  A  and  S2l  for 
fixture  B  are  approximately  the  same,  one  can  choose  either  of  the  phase  angles.  Some 
confusion  may  arise,  if  the  scattering  parameters  of  both  fixture  A  and  B  are  specified 
independently  of  the  final  device  (51  ohm  resistor)  parameters.  In  this  case,  the  ap¬ 
proximate  phase  shift  through  otic  fixture  can  be  computed,  or  determined  exper- 


imentally,  and  used  to  choose  the  proper  phase  of  S2I  for  fixture  A  and  S21  for  fixture 

B. 


In  Chapter  V,  the  de-embedding  procedure  with  symmetrical  fixtures  was  carried 
out,  using  as  the  device  a  commercial  bandpass  filter  (8-12  GHz).  Standard  hardware 
couplers  were  added  to  represent  the  fixtures.  Since  the  fixtures  were  symmetrical,  the 
de-embedding  procedure  was  relatively  straightforward.  As  may  be  seen  in  Figs.  26  to 
2S.  the  de-embedded  scattering  parameters  of  the  filter  were  very  close  to  the  expectation 
values,  i.e.  of  the  scattering  parameters  of  filter  without  fixtures.  It  may  also  be  noted 
that  the  measured  S-parameters  of  the  embedded  filter  (fixtures  plus  filter)  were  very 
similar  to  those  of  the  filter  alone,  except  of  course  for  phase;  indicative  of  well-matched 
coupling  hardware.  Therefore,  in  this  task,  the  de-embedded  results  were  found  to  be  in 
good  agreement  with  the  expected  values,  as  shown  in  Figs.  26  through  28. 
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APPENDIX  A.  THE  OUTPUTS  OF  THE  NETWORK  ANALYZER  AND 
THEIR  REFLECTION  COEFFICIENTS  FOR  FIXTURE  A  AND  B 


T 

C 

C\J 


Figure  33.  Measurement  of  the  Load  1  (Fixture  A) 


Table  5.  THE  MEASUREMENT  OF  THE  LOAD  KFIXTURE  A) 


Freq.(  MHz) 

<JB  (loss) 

Angle!  dee.) 

r. 

2000.0 

3.07 

77  7 

0.65  +  jo. 27 

2150.0 

3.82 

-63.1 

0.29  -  jo. 57 

2300.0 

3.70 

-153.9 

-0.58  -  jo. 28 

2450.0 

3.46 

126.4 

■■desssbh 

2600.0 

3.20 

60.5 

HHHEERflBDH 

2750.0 

3.33 

-1.7 

■■DESBOSai 

2900.0 

4.04 

-69.9 

0.22  -  jO. 59 

3050.0 

5.39 

-152.S 

-0.4S  -  jo. 25 

3200.0 

6.17 

107.6 

-0.15  + jo.47 

3350.0 

5.32 

17.1 

HHEE33DZ9I 

3500.0 

4.65 

-51.0 

■■BSEBRESH 

3650.0 

3.97 

-111.9 

■■bbhqsb 

3800.0 

4.89 

-175.9 

HHSEflRIB9i 

3950.0 

7.73 

104.5 

■HBSBtEH 

4100.0 

7.36 

-17.1 

0.41  -  jo. 13 

4250.0 

6.70 

-121.4 

-o.24  -  jo. 39 

4400.0 

4.59 

170.4 

mmseowm 

4550.0 

4.12 

114.2 

-0.26+  jo.57 

4700.0 

4.88 

58.9 

j^HKESDIEEK 

4850.0 

7.91 

-9.1 

0.4o  -  jo. 06 

5000.0 

8.42 

-137.9 

-0.28  -  jO.25 

57 


Dfinector  !  *r  ! 


Figure  34.  Measurement  of  the  Load  2  (Fixture  A) 


Table  6.  THE  MEASUREMENT  OF  THE  LOAD  2  (FIXTURE  A> 


Freq.(MHz) 


2000.0 


150.0 


2300.0 


2450.0 


2600.0 


2750.0 


2900.0 


3050.0 


3200.0 


3350.0 


3500.0 


3650  0 


3800.0 


3950.0 


4100.0 


4250.0 


4400.0 


4550.0 


470(3.0 


4850.0 


5000.0 


dB  (loss* 


Anglet  deg.) 


9S.2 


29.5 


-54.2 


-119.0 


o.l  1  +  in.  76 


0.62  ■+•  jo.3  5 


0.42  -  jo.49 


0.67  -t-  j().06 


0.39  -  jO.54 


-0.4S-t-j0.lS 


0.23  + jO.42 


0.01  -  jO. 5 2 


-0.41  -  jo.42 


-0.63  +  jo.07 


-0.27  f  jO. 53 


0.31+ jO.  34 


o. IS  -  jO.33 


.0,42  -  jo. i >4 


-0.24*  jo.  5 1 


0.30  +  j0.56 


; 

toad  3  > 


Figure  35.  Measurement  of  the  Load  3  (Fixture  A) 


Table  7.  THE  MEASUREMENT  OF  THE  LOAD  3  (FIXTURE  A) 


Prcq.(MHz)  d B  floss)  Anciefdeg.) 


.14 


2.21 


2()()u.o 


2150.0 


2450.0 


2600.0 


2750.0 


2000. 0 


3200.0 

3.89 

3350.0 

3.89 

3500.0 

4.15 

3650.0 

5.02 

3800.0 

8.22 

3950.0 

8.95 

4100.0 

6.37 

4250.0 

4.35 

4400.0 

3.49 

4550.0 

3.37 

4700.0 

3.90 

4850.0 

5.55 

5000.0 

6.59 

-66.5 


-120.7 


175.9 


99.0 


9.5 


-79.5 


-147.2 


157.8 


1 08.4 


58. 


1.4 


-82.6 


-0.78  +  jO.O 

/ 

-0.32  - 

jo.  7 

1 

0.42 

jo.6 

1 

_ 

0.69  . 

jo.  06 

0.1"  - 

jo. 6 

■A 

-0.55  - 

wm 

-o.3<»- 

i'».5 

0.43  - 

jit. 46 

0.62  - 

ill.  lb 

0.2<  - 

jo.5 

6 

-0.32  - 

j0-5 

3 

-0.56  + j0.04 

-0.06  +  jO.  3  8 

0.35  +  j0.06 

0.09  - 

jO.47 

-0.62  — 

j0-2 

5 

-0.2 1  +  jO. 64 

0.33  -i- 

jo. 55 

0.53  + 

j(l.() 

i 

0.06  - 

jO. 46 

Figure  37.  Measurement  of  the  Load  2  (Fixture  B) 


Table  9.  THE  MEASUREMENT  OF  THE  LOAD  2  (FIXTURE  B) 


Freq.(  MHzi 

dii  (loss) 

Ancle!  dec.) 

r,„ 

2000.0 

2.32 

-26.9 

0.68  -  jo. 35 

2150.0 

2.46 

-91.9 

-0.02  -  jo.75 

230o. o 

2.44 

-155.5 

-0.69  -  jo. 31 

245o.  o 

2.51 

147.3 

-0.63  —  j0.4i  i 

2600.0 

2.41 

96.6 

-0.09  + jo.75 

2750.0 

2.25 

50. 1 

0.50  +  jo.59 

2900.O 

2.35 

6.9 

0.76  +  jo.  09 

3050.0 

2.2S 

-35.9 

0.62  -  jO.45 

320o.ii 

2.51 

-78. 5 

0.15  -  i0.73 

335o.o 

2.69 

-123.9 

-0.41  -  j0.6! 

3500.0 

3.25 

-175.1 

-0.69  -  jO.06 

3650.0 

4.27 

126.2 

-0.36  + jO. 49 

38oo.O 

5.24 

56.1 

0.31  +  jO.45 

3950.0 

5.38 

-20.6 

0.50  -  j0. 19 

4100.0 

4.49 

-90.2 

o.oo  -  jO. (4  • 

4250.0 

3.49 

-147.9 

-0.57  -  j0.3(' 

4400.0 

2.76 

163.0 

-0.7o  — ]0.21 

4550.0 

2.43 

I  IS. 9 

-0.37 -j.  (.of  i 

4700.0 

2.36 

77.5 

0.16+ jO. 74 

4850.0 

2.55 

36.2 

0.60  +  jO. 44 

5000.0 

3.32 

-8.6 

0.67  -  jO.10 

Figure  38.  Measurement  of  the  Load  3  (Fixture  B) 


Table  10. 


THE  !\ 


EASUREMENT  OF  THE  LOAD  3  (FIXTURE  B) 


Freq.(  MHz) 

dB  (loss) 

r. 

2000.0 

1.41 

7S.5 

■HBBnraa 

2150.0 

1.71 

2S.1 

( ». “2?  -  id.39 

2300.0 

2.16 

-2o.2 

0“O  -  jt  >.34 

2450.0 

2.4S 

-84.3 

o.o"  -  jo.75 

2600.0 

2.69 

-145.9 

-0.61  -  ]0.4l 

275o.o 

r  j 

be 

*s>l 

152.4 

-0.64-  jO. 3.3 

2900.0 

2. 85 

94.8 

-0.06  - 10.72 

3(150.0 

2.75 

43.2 

0.53  ]0.5o 

320O.0 

2.77 

-2.9 

o.7*  -  jo.o4 

3350.0 

2.68 

-44.3 

0.53  -  jo.5 1 

3500.0 

2.6S 

-S2.4 

o.lu  -  io.73 

3650.0 

1.93 

-121.1 

-0.41  -  jo.69 

3800.0 

1.91 

-162.1 

-0.76  -  io.25 

3950.0 

2.02 

157.0 

-0.73  ■+■  jo. 3 1 

4100.0 

2.42 

113.2 

-0.3O jo.7o 

4250.0 

3.21 

64.6 

o.  3o  +  jo.  62 

4400.0 

4.49 

7.5 

0.59  |0.0S 

4550.0 

5.70 

-60.8 

0.25  -  jo. 45 

4700.0 

5. SO 

- 1 34.9 

4850.0 

4. S') 

158.8 

mssAzsH 

5000.0 

3.72  105.1 

-U.17-1- ju.63 

1 


APPENDIX  B.  THE  S-PARAMETERS  OF  THE  FIXTURES 


Table  11.  THE  S  -  PARAMETERS  OF  FIXTURE  A 


Freq.(MHz) 

^12 

A  3 

mmm 

-0.130  +  j0.0S2 

O.S37  -  jO.  143 

0.132  4  j0.029 

150.0 

-0.064  4  jO.  153 

0.668  -  jO. 499 

0.139  -  j0.019 

2300.0 

0.026  4  jO.  165 

0.355  -  jO.738 

0.130  -  jO.OSO 

“> 

450.0 

0.099  4  jO.  125 

0.021  4  j0.802 

0.100  -  jO.  147 

600.0 

0.141  4-  j0.069 

0.380  4  jO.701 

0.034  -  jO.  17S 

2750.0 

0.151  -  j0.007 

0.653  +  j0.443 

-0.018  -  j0.180 

2900.0 

0.091  -  j0.063 

0.777  4  j0.0S7 

-0.032  -  jO.  128 

3050.0 

0.048  -  jO.055 

0.717  -  j0.287 

-0.003  -  j0.074 

3200.0 

0.042  -  j0.026 

0.485  -  j0.576 

0.060  -  j0.037 

3350.0 

0.098  -  j0.025 

0.151  -  j0.720 

0.180  -  j0.031 

J 

500.0 

0.137  -  jO.l  13 

0.199  4  j0.679 

0.2S8  -  jO.OSS 

3650.0 

0.054  -  j0.265 

0.470  4  j0.482 

0.330  -  jO.  184 

3S00.0 

-0.148  -  j0.261 

0.579  4  j0.2 14 

0.327  -  jO.254 

3950.0 

-0.258  -  jO.  163 

0.625  -  j0.036 

0.321  -  jO.239 

4100.0 

-0.334  4  J0.061 

0.581  -  j0.305 

4250.0 

-0.256  4  jO.288 

0.318  -  j0.506 

0.442  -  jO.  181 

4400.0 

0.174  4  jO.317 

0.008  -  j0.477 

0.713  -  j0.044 

4550.0 

-0.502  -  j0.872 

0.225  -  j  1 .580 

1.743  -  j  1 .374 

4700.0 

-0.064  -  j0.53 1 

0.670  -  jO.309 

0.034  -  j0.761 

4850.0 

0.183  4  j0.229 

0.651  4  j0.057 

0.137  -  j0.304 

Table  12.  THE  S  -  PARAMETERS  OF  FIXTURE  B 


Frcq.(MHz) 

Du 

Du 

Dr. 

| 

-0.060  +  j0. 107 

0.342  -  jO. 834 

0.036  -  j0.076 

2150.0 

-0.007  +  jO.  1 16 

0.O0S  -  jO.OSl 

2300.0 

0.027  4-  jO.  100 

0.230  +  j0.850 

-0.032  -  jO.085 

2450.0 

0.049  +  j0.091 

0.475  +  jO.  723 

-0.083  -  j.082 

2600.0 

0.0S0  +  j0.084 

0.669  +  jO.533 

-0.115  -  jO.075 

2750.0 

0.124  +  j0.063 

0.792  -  jO. 29 1 

-0.125  -  jO.062 

2900.0 

0.152  +  j0.007 

0.831  +  jO.028 

-0.136  -  j0.029 

3050.0 

0.158  -  jO.073 

0.788  -  jO. 229 

-0.125  -  j0.007 

3200.0 

0.106  -  jO.  143 

0.669  -  j0.451 

-0.110  +  jO.016 

3350.0 

0.026  -  jO.168 

0.494  -  jO.630 

-0.064  +  j0.036 

3500.0 

-0.035  -  jO.  149 

0.268  -  jO.741 

0.016  +  j0.059 

3650.0 

-0.092  -  jO.  112 

0.033  -  j0.790 

0.108  +  jO.019 

3800.0 

-0.076  -  j0.072 

0.230  4-  j0.762 

0.165  +  jO.018 

3950.0 

-0.069  -  j0.097 

0.468  4-  jO.633 

0.261  4-  jO.019 

4100.0 

-0.137  -  jO.  157 

0.636  +  j0.412 

0.389  4-  jO.0 1 4 

4250.0 

-0.331  -  jO.  130 

0.681  +  jO.  135 

0.543  4-  j0.018 

4400.0 

-0.572  4-  jO.  161 

0.478  -  jO.  116 

0.808  4-  j0.024 

4550.0 

-0.114  +  j  1 . 1 7 1 

1.062  4-  jO.589 

1.625  -  jO.527 

4700.0 

0.230  -  jO.  171 

0.577  -  j0.759 

0.163  -  j0.968 

4850.0 

-0.065  +  jO.  166 

0.306  -  j0.676 

0.219  -  j0.546 

5000.0 

0.124  +  jO.326 

0.091  -  j0.685 

0.291  -  j0.446 

APPENDIX  C.  THE  S-PARAMETERS  OF  THE  EMBEDDED  AND 
DE-EMBEDDED  51  OHM  RESISTOR 


Table  13.  THE  THEORETICAL  S  -  PARAMETERS  OF  THE  51 
OHM  RESISTOR  IN  DB  UNITS. 


Freq.(MHz) 

s„ 

S:2 

2000.0 

6.11 

4.27 

4.27 

6.11 

2150.0 

5.82 

4.37 

4.37 

5.82 

2300.0 

5.54 

4.47 

4.47 

5.54 

2450.0 

5.28 

4.58 

4.58 

5.28 

2600.0 

5.02 

4.69 

4.69 

5.02 

2750.0 

4.78 

4.81 

4.81 

4.78 

2900.0 

4.55 

4.93 

4.93 

4.55 

3050.0 

4.34 

5.05 

5.05 

4.34 

3200.0 

4.13 

5.17 

5.17 

4.13 

3350.0 

3.94 

5.30 

5.30 

3.94 

3500.0 

3.76 

5.43 

5.43 

3.76 

O' 

o 

b 

3.58 

5.56 

5.56 

3.58 

3S00.0 

3.42 

5.69 

5.69 

3.42 

3950.0 

3.27 

5.82 

5.82 

3.27 

4100.0 

3.12 

5.95 

5.95 

3.12 

4250.0 

2.99 

6.09 

6.09 

2.99 

4400.0 

2.86 

6.22 

6.22 

2.86 

4550.0 

2.73 

6.36 

6.36 

2.73 

4700.0 

2.62 

6.49 

6.49 

2.62 

4S50.0 

2.51 

6.03 

6.63 

2.51 

5000.0 

2.41 

0.-6 

6.76 

2.41 

Table  14.  THE  S  -  PARAMETERS  OF  BOTH  FIXTURE  AND  51 
OHM  RESISTOR 


Freq.(MHz) 

A/u 

Mn 

.U:, 

Mu 

■ 

11.14 

4.40 

4.60 

6.45 

2150.0 

11.3S 

4.64 

4.71 

6.27 

2300.0 

10.18 

4.90 

4.80 

6.05 

2450.0 

8.47 

5.36 

5.35 

5.86 

2600.0 

7.13 

5.87 

5.84 

5.51 

2750.0 

6.35 

6.30 

6.25 

4.99 

2900.0 

6.21 

6.99 

6.99 

4.68 

3050.0 

6.51 

7.54 

7.49 

4.48 

3200.0 

7.33 

7.32 

7.25 

4.43 

3350.0 

8.92 

7.04 

6.65 

4.94 

3500.0 

11.29 

7.10 

5.81 

5.75 

3650.0 

12.01 

6.46 

6.31 

6.18 

3800.0 

11.59 

7.30 

7.25 

6.91 

3950.0 

9.15 

6.45 

6.90 

7.83 

4100.0 

6.27 

6.35 

6.56 

9.07 

4250.0 

5.53 

6.00 

6.25 

1.78 

4400.0 

6.07 

5.61 

5.85 

12.72 

4550.0 

7.54 

5.30 

5.50 

12.53 

4700.0 

9.91 

5.34 

5.39 

10.00 

4850.0 

10.46 

5.70 

5.66 

7.82 

5000.0 

11.95 

6.10 

6.05 

6.53 

Table  15.  THE  DE  -  EMBEDDED  OUTPUT  OF  51  OHM  RE 
SISTOR  DEVICE  UNDER  TEST. 


Freq.(MHz) 

s„ 

sl2 

6.35 

3.62 

3.82 

2150.0 

5.72 

3.97 

3.64 

2300.0 

5.12 

3.97 

3.97 

2450.0 

4.57 

3.18 

3.17 

2600.0 

4.22 

3.39 

3.36 

2750.0 

3.61 

3.66 

3.61 

2900.0 

3.87 

4.27 

4.27 

3050.0 

3.64 

4.91 

4.86 

3200.0 

3.19 

4.51 

4.44 

3350.0 

3.03 

4.16 

3.77 

3500.0 

3.16 

4.72 

4.43 

3650.0 

3.62 

4.65 

4.50 

3S00.0 

3.91 

4.56 

4.51 

3950.0 

3.35 

5.38 

4.83 

4100.0 

2.66 

5.06 

5.27 

4250.0 

3.45 

5.71 

5.46 

4400.0 

4.11 

5.83 

6.07 

4550.0 

3.17 

6.36 

6.56 

4700.0 

2.25 

6.73 

6.68 

4850.0 

2.08 

6.73 

6.69 

5000.0 

2.11 

6.96 

6.91 

APPENDIX  D.  THE  S-PARAMETERS  OF  THE  FIXTURE,  FILTER 
EMBEDDED,  FILTER  WITHOUT  FIXTURES  AND  FILTER 

DE- EM  BEDDED 


Table  16.  THE  S  -  PARAMETERS  OF  THE  FIXTURE 


Freq.(MHz) 

s„ 

S,i 

$22 

6000.0 

-0.08  +  jO.  1 7 

0.81  +  j0.55 

-0.09  +  jO.  1 3 

6320.0 

-0.09  +  j0.04 

0.76  -  j0.6 1 

-0.17  +  jo.  10 

6640.0 

-0.06  4-  j0.04 

-0.39  -  j0.90 

-0.03  +  jO.  1 2 

6960.0 

-0.05  4-  j0.03 

-0.98  4-  j0.16 

-0.08  4-  j0.07 

72S0.0 

-0.19  +  j0.03 

-0.08  4-  j0.97 

-0.1S  +  jO.  15 

7600.0 

-0.13  +  jO.  14 

0.92  4-  jO.36 

-0.05  +  jO.  18 

7920.0 

0.01  4-  j0.07 

0.57  -  j0.81 

0.01  +  j0.ll 

8240.0 

-0.07  +  jO.OO 

-0.63  -  j0.76 

0.U2  4-  j0.ll 

S560.0 

-0.12  4-  j0.04 

-0.90  4-  j0.42 

0.03  -  jO.Ol 

88S0.0 

-0.10  4-  j0.07 

0.18  4-  j0.97 

-0.11  -  j0.03 

9200.0 

0.00  +  j0.19 

0.96  4-  j0.06 

-0.03  4-  jO.  20 

9520.0 

0.10  4-  j0.14 

0.33  -  jO.9 1 

0.11  4-  jO.ll 

9840.0 

0.11  -  j0.04 

-0.81  -  j0.55 

-0.02  -  j0.09 

10160.0 

0.12  -  jO.  19 

-0.74  4-  j0.66 

-0.03  -  j0.03 

10480.0 

-0.03  -  jO.  13 

0.42  4-  jO.88 

0.03  4-  j0.02 

10800.0 

-0.02  -  jO.  14 

0.96  -  jO.  18 

-0.02  -  j0.09 

11120.0 

-0.18  -  j0.l2 

0.06  -  30.98 

-0.16  4-  j0.03 

11440.0 

-0.14  -  j0.03 

-0.91  -  jO.33 

-0.08  4-  j0.02 

11760.0 

-0.02  -  jO.  14 

-0.57  4-  j0.79 

-0.04  -  j0.04 

12080.0 

-0.04  -  j0.04 

0.63  4-  j0.75 

0.04  4-  j0.06 

12400.0 

-0.02  -  jO.10 

0.85  -  j0.45 

-0.02  -  j0.04 
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Table  17.  THE  S  -  PARAMETERS  OF  FILTER  +  FIXTURE 

Freq.(MHz)  Sn  S12 


6000.0 

-0.79  -  j0.45 

0.00  +  jO.OO 

0.54  +  jO.75 

6320.0 

0.S9  -  jO.  14 

0.00  +  jO.OO 

-0.36  -  jO.S2 

6640.0 

-0.65  -t-  j0.62 

0.01  +  jO.Ol 

-0.02  +  j  0 . 9  0 

6960.0 

0.11  -  j0.9 1 

0.01  +  jO.Ol 

0.05  -  jO.S9 

72S0.0 

0.62  +  j0.64 

0.01  +  jO.Ol 

-0.24  +-  j0.90 

7600.0 

-0.S0  +  j0.45 

-0  02  +  j0.02 

0.02  -  j0.91 

7920.0 

0.11  -  j0.24 

-0.30  +  j0.86 

0.00  +  j0.25 

8240.0 

-0.13  4-  j0.20 

-0.34  +  j0.86 

-0.01  -  jO.  1 5 

8560.0 

-0.01  +-  j0.36 

-0.82  +  j0.24 

0.15  +  j0.27 

SSSO.O 

0.01  +  j0.04 

-0.61  -  j0.69 

-0.0 1  -  j0.05 

9200.0 

0.31  +  jO.Ol 

0.56  -  j0.65 

0.03  +  j0.31 

9520.0 

0.27  +  j0.29 

0.76  +  jO.41 

-0.28  -  j0.06 

9S40.0 

-0.16  +  jO.22 

-0.48  +  j0.76 

0.07  -  jO.23 

10160.0 

-0.06  +  jO.23 

-0.73  -  j0.50 

0.00  +  j0.29 

10480.0 

-0.06  -  j0.26 

0.47  -  j0.76 

-0.26  -  j0.03 

10800.0 

0.21  -  j0.02 

0.79  +  j0.46 

-0.05  -  jO.  19 

11120.0 

-0.17  -  j0.07 

-0.45  +  jO.78 

-0.14  -  j0.02 

11440.0 

-0.14  +  jO.OO 

-0.83  -  jO.37 

0.14  +  jO.Ol 

11760.0 

0.19  +  j0.07 

0.16  -  j0.86 

0.17  +  j0.05 

12080.0 

0.25  -  j0.24 

0.72  -  j0.28 

-0.48  -  j0.13 

12400.0 

-0.15  +  j0.70 

0.74  +  jO.18 

0.59  +  jO.  12 

* 


Table  18.  THE  S  -  PARAMETERS  OF  FILTER  WITHOUT  FIX 
TURES  IN  DB  UNITS 


Frcq.(MHz) 

s.. 

sa 

sn 

^22 

6000.0 

1.57 

83.91 

88.79 

0.62 

6320.0 

1.96 

91.38 

105.37 

1.69 

6640.0 

0.47 

75.81 

7S.12 

0.19 

6960.0 

0.87 

69.94 

69.39 

0.85 

7280.0 

0.77 

54.40 

51.94 

0.86 

7600.0 

0.70 

29.99 

31.60 

1.12 

7920.0 

10.17 

0.57 

0.50 

9.15 

S  240.0 

11.32 

0.15 

0.16 

12.00 

8560.0 

8.62 

0.45 

0.45 

10.05 

S8S0.0 

15.77 

0.20 

0.08 

13.81 

9200.0 

9.66 

0.40 

0.20 

9.88 

9520.0 

10.97 

0.25 

0.33 

11.81 

9840.0 

9.65 

0.27 

0.66 

9.11 

10160.0 

8.39 

0.45 

0.40 

10.81 

10480.0 

12.78 

0.15 

0.20 

12.01 

10800.0 

11.35 

0.43 

0.45 

11.87 

11120.0 

13.64 

0.30 

0.20 

14.39 

11440.0 

10.28 

0.31 

0.46 

16.09 

11760.0 

15.44 

0.60 

0.41 

15.84 

12080.0 

7.35 

1.40 

1.10 

9.98 

12400.0 

8.85 

1.45 

0.96 

11.51 

Table  19.  THE  DE-EMBEDDED  OUTPUT  OF  FILTER  UNDER 
TEST 


m 

M,V»| 

-  M 
■sn 

wi 

wl 

*M*1 

•Ml 


•M2I 

yi 


M 


>*  I*1 

v* iii 

•  vi 
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m 

0.61 


Freq.(MHz) 

5,. 

S,2 

■Sji 

6000.0 

1.58 

37.03 

90.64 

0.28 

6320.0 

2.03 

41.18 

91.12 

1.15 

6640.0 

1.11 

41.62 

76.82 

0.14 

6960.0 

2.32 

59.94 

71.12 

0.90 

72S0.0 

0.68 

51.99 

53.34 

1.21 

7600.0 

0.21 

30.06 

31.45 

0.92 

7920.0 

10.50 

1.27 

1.51 

9.64 

S  240.0 

10.99 

0.59 

0.60 

13.25 

8560.0 

7.97 

1.25 

1.18 

10.17 

SSSO.O 

12. S3 

1.12 

0.55 

10.50 

9200.0 

9.57 

1.36 

0.70 

10.77 

9520.0 

8.80 

0.35 

0.94 

9.70 

9S40.0 

10.83 

0.11 

0.81 

8.38 

10160.0 

7.58 

0.74 

0.90 

10.05 

104S0.0 

11.96 

0.62 

0.75 

9.76 

10S00.0 

12.39 

0.67 

0.73 

9.99 

11120.0 

11.49 

0.14 

0.20 

9.78 

1 1440.0 

8.94 

0.3S 

0.76 

15.48 

11760.0 

13.73 

0.76 

0.94 

14.58 

12080.0 

8.20 

2.11 

1.41 

6.85 

12400.0 

8.41 

1.77 

1.38 

7.79 

wmmmmm 


APPENDIX  E.  FORTRAN  PROGRAM  FOR  DE-EMBEDDING  THE  51 

OHM  RESISTOR 


DIMENSION  X(200) ,  MDB 11(200) ,PSE11(200) 

MDB12( 200), PSE 12(200) 

MDB21(200) ,PSE21(200) 

MDB22( 200), PSE 22(200) ,MDB222(200) 

REAL  MDB 1 1 , MDB 1 2 , MDB  2 1 , MD  B  2  2 , MDB  2  2  2 
MAG 1 1 , MAG 12 , MAG 2 1 , MAG22 
IM11, IM12 , IM21 , IM22 
COMPLEX  S11,S12,S21,S22 ,A12A21 , B 12B21 

All ,A12 , A21 , A22 , AAI21 , A12B12 ,A21B21 
B11,B12,B21,B22,BB121 
Mil ,M12 ,M21 ,M22 
E,FF,G,H,N,AABB 
PI=4. *ATAN( 1.  ) 

DO  10  1=1,21 

F=18.  5*10. **8+1. 5*10. **S*rL0AT( I) 

X( I)=F/ 1000000. 

CALL  FIX1( All ,A12,A21 ,A22,A12A21 , AMAG12 ,PSEA12 ,F) 
CALL  TMFS(M11 ,M12 ,M21 ,M22,F) 

CALL  FIX2(B11 ,B12 ,B21 ,B22 ,B12B21 ,BMAG12 ,PSEB12 ,F) 


TMAG=AMAG 1 2*BMAG 1 2 
TPSE=PSEA12+PSEB12 
XRE=TMAG*COS( TPSE ) 

XIM=TMAG*SIN(TPSE) 

A12B 12=CMPLX( XRE , XIM) 

A21B21=A12B12 

E=M12/(A12B12) 

FF=M21/( A21B21) 

G=(M11-A11)/A12A21 

H=(M22-B22)/B12B21 

N=( 1. +A22*G)*( 1.  +B1 1*H) -A22*B 1 1*E*FF 


S 1 1=( G*( 1 .  +B 1 1*H ) - B 1 1*E*FF ) / N 

S12=E/N 

S21=FF/N 

S22=(H*( 1. +A22*G)-A22*E*FF)/N 


QIM1=AIMAG( Sll) 


QRE1=REAL(S11) 

QIM2=AIMAG( S12) 

QRE2=REAL(S12) 

QIM3=AIMAG( S21) 

QRE3=REAL(S21) 

QIM4=AIMAG( S22) 

QRE4=REAL(S22) 

FIN'D  A  MAGNITUDE  AND  PHASE( ANGLE) 

MAG  1 1=C ABS ( S 1 1 ) 

MAG12=CABS(S12) 

MAG2 1=CABS( S2 1 ) 

MAG22=CABS( S22) 

MDB11( I)=ABS( 20*ALOG10( MAGI 1) ) 

MDB12( I )=ABS( 20*ALOG10( MAG12) ) 

MDB21( I)=ABS(20*ALOG10(MAG21) ) 

MDB22( I)=ABS( 20*ALOG10( MAG22 ) ) 
PSE11(I)=ATAN2(QIM1,QRE1)*180. /PI 
PSE12( I )=ATAN2( QIM2 ,QRE2)*180. /PI 
PSE21(I)=ATAN2(QIM3,QRE3)*180. /PI 
PSE22(  I)=ATAN2(QIM4,QRE4)*180. /PI 

WRITE(6, 101)  X( I) ,MDB11( I) ,MDB12( I) ,MDB21( I) ,MDB22( I) , 
1  PSE11C I) ,PSE12( I) ,PSE21( I) ,PSE22( I) 

1  F0RMATCF6.  0,8(1X,F7.  2)) 

CONTINUE 

END 


SUBROUTINE  TMFS  (Dll ,D12,D21,D22,F) 

VoWr*')V****')’c******iWfiV***')V*iWf***iWr>ViV7'r,iV 


COMPLEX  011,012,021,022 
PI=4. *ATAN( 1.  ) 

INPUT  DATA  FILE  IS  DISK  9 
(  BOTH  FIXTURE  +  DEVICE  ) 

READ( 9,21)DB11 , ANG11 ,DB21 , ANG21 ,DB12 , ANG12 ,DB22 , ANG22 
FORMAT( 8(F7. 2)) 

Vll=10. **( -DB11/20. ) 

V12=10.  **( -DB12/20.  ) 

V21=10.  **( -DB21/20.  ) 

V22=10.  **( -DB22/20.  ) 

RAD11=PI/ 180.  *ANG11 
RAD12=PI/180.*ANG12 
RAD21=PI/180.  *ANG21 
RAD22=PI/180.*ANG22 

QRE 1 1=V 1 1*C0S( RAD  11) 

QRE12=V12*COS(RAD12) 

QRE21=V21*C0S(RAD21) 


D 
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in 

|V 
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v. 
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QRE22=V22*C0S( RAD22 ) 
QIM1 1=V 1 1*S IN( RAD 11) 
QIM12=V12*SIN( RAD12) 
QIM21=V21*SIN(RAD21) 
QIM22=V22*SIN( RAD22) 


D11=CMPLX(QRE11,QIM11) 
D12=CMPLX(QRE12 ,QIM12) 
D2 1=CMPLX( QRE2 1 , QIM2 1 ) 
D22=CMPLX( QRE22 ,QIM22) 


RETURN 

END 


SUBROUTINE  FIX1( Sll , S12 , S21 , S22 ,SS , SMAG12 , SPSE12 ,F) 

iV  ~'c  V?  Vc  Vc  V?  it  ?V  V?  it  -jV  Vc  V?  Vc  :  V  Vc  Vc  Vc  Vc  Vc  Vr  Vr  Vc  V?  Vc  -/V  -aWr  -,V  ?V  Vr  iV  Vc  -.V  Vc-  Vc  *  Vo V  *  Vc*  -  V  Vofc  Vr  Yr  Vr  Vr  Vo Wf  Yr 


FOR  FIXTURE  A 


COMPLEX  A(3,3),H(3,3),B(1,3),C(3),D(3,3),SS121 
COMPLEX  311,812,321,322,33 


CALL  DATA1(A,B,F) 


CALL  CPX(A) 


DO  10  1=1,3 

C( I)=A( 1,I)*B( 1, 1)+A(2,I)*B( 1,2)+A(3,I)*B( 1,3) 
CONTINUE 


S11=C(1) 

S22=C( 2) 

S  S 1 2 1=C ( 3) 
SS=S11*S22-SS121 


SIM=AIMAG( SS) 
SRE=REAL( SS) 
SMAGSS=CABS(SS) 
PSESS=ATAN2( SIM, SRE) 
SPSE12=0.  5*PSESS 
SPSE21=SPSE12 
SMAG12=SMAGSS**0. 5 
SMAG21=SMAG12 


S12=(S11*S22-SS121)**0. 5 

S21=S12 

RETURN 

END 


*********************** 


SUBROUTINE  DATA1(D,B,F) 
*********************** 


COMPLEX  D(3,3) ,B( 1,3) ,RIN1 ,RIN2 ,RIN3 ,RL1 ,RL2 ,RL3 


CALL  RLOAD( RL1 , RL2 , RL3 , F ) 

CALL  RIN(RIN1,RIN2,RIN3) 

D( 1 , 1)=1. /RL1 
D( 1 , 2)=1. /RL2 
D(  1 , 3)  =  1.  /RL3 
D( 2 , 1)=RIN1 
D( 2 , 2)=RIN2 
D( 2 , 3)=RIN3 
D(  3 , 1  )=(  - 1.  ,.0) 

D(  3 , 2)=( - 1.  0) 

D(  3 , 3 )=( - 1.  , .  0 ) 

B( 1 , 1)=RIN1/RL1 
B( 1 , 2 )=RIN2/RL2 
B( 1 ,3)=RIN3/RL3 
RETURN 
END 


is  is  is  is  is  isisisisisis  is  is  isisisisisisisisisisisisisisis  is  is 

SUBROUTINE  RIN(RIN1 ,RIN2 ,RIN3) 

isisisisisisisisisicisisisitisisisisisisisisisisisisisisisis 


COMPLEX  RIN1 ,RIN2 ,RIN3 

CALL  F1DAT(DB111 , ANG111 ,DB112 , ANG112 ,DB113 , ANG113) 

PI=4. *ATAN( 1.  ) 

INPUT  DATA  FILE  IS  DISK8 
(FIXTURE  1  ONLY) 

Vlll=10.**(-DBlll/20.  ) 

VI 12=10.  **( -DB 112/20.  ) 

VI 13=10. **( -DB113/20.  ) 

RAD111=PI/180.  *ANG111 
RADI 12=PI/ 180.  *ANG112 
RAD113=PI/180.  *ANG113 

QRE111=V111*C0S(RAD111) 

QRE112=V112*COS(RAD112) 

QRE113=V113*C0S(RAD113) 

QIM111=V111*SIN(RAD111) 

QIM112=V112*SIN(RAD1 12) 

QIM113=V113*SIN(RAD113) 

RIN1=CMPLX(QRE111 , QIM111) 

RIN2=CMPLX(QRE112 ,QIM112) 

RIN3=CMPLX(QRE113,QIM113) 

RETURN 

END 


******************************************************** 
SUBROUTINE  F1DAT(DB111 , ANG1 1 1 , DB 1 12 , ANG112 ,DB113 ,ANG113) 

*  *  *  *  *  isisisis  irisisiciSisisisisis  isisisis  is  *  V:  is  is  *  *  is  is  is  is  is  is  i:  iris  it  isisisisis  isisis  is  is  is  is  is  is  is 


INPUT  DATA  FILE  IS  DISK  8 
(For  FIXTURE  A) 

READ( 8 , 20)DB11 1 , ANG111 ,DB112 , ANG112 ,DB113 , ANG113 
FORMAT( 6(F7. 2) ) 

RETURN 

END 


**Vf*-*Vf’*«A,*V?*****<;WrVf*Vr***«>Vv-*Vf*<;V*.Wf*Vf***Vr**y?V?Vc,;W?***Vc 

SUBROUTINE  FIX2( Sll , S12 , S2I , S22 , SSS , SMAG12 , SPSE12 ,F) 

*****VrVc**V{,Vf****yr>Wf*Vr*V,-,;,f**V?';V******yrVr,5W»'*Vr*Vf*yrjVVfVf** 


FOR  FIXTURE  2 


COMPLEX  A(3,3),H(3,3),B(1,3),C(3),D(3,3),SS121 
COMPLEX  311,812,321,322,333 

CALL  DATA2( A, B ,F) 

CALL  CPX(A) 

DO  10  1=1,3 

C(I)=A(1,I)*B(1,1)+A(2,I)*B(1,2)+A(3,I)*B(1,3) 

CONTINUE 

S11=C(1) 

S22=C( 2) 

SS12l=C( 3) 

SSS=S11*S22-SS121 


SIM=AIMAG(SSS) 
SRE=REAL( SSS) 
SMAGSS=CABS( SSS) 
PSESS=ATAN2( SIM.SRE) 
SPSE12=0. 5*PSESS 
SPSE21=SPSE12 
SMAG12=SMAGSS**0. 5 
SMAG21=SMAG12 


S12=(S11*S22-SS121)**0.  5 

S21=S12 

RETURN 

END 


SUBROUTINE  DATA2(D,B,F) 

*Vt***'********TWr*'i,r**Vf'>lfVfVf 


COMPLEX  D(3,3) ,B( 1 ,3) ,RIN1 ,RIN2 ,RIN3 ,RL1 ,RL2 ,RL3 
CALL  RL0AD(RL1 ,RL2 ,RL3 ,F) 

CALL  R2IN(RIN1 ,RIN2 ,RIN3) 


D(  1 , 1)  =  1.  /RL1 
D( 1 , 2)=1. /RL2 
D( 1 , 3)=1. /RL3 
D( 2 , 1)=RIN1 
D( 2 , 2)=RIN2 
D( 2 , 3)=RIN3 
D(  3 , 1)=( - 1.  0) 

D( 3 , 2)=( -1. , . 0) 
D(  3 , 3 )=(  - 1.  ,.0) 
B( 1 , 1)=RIN1/RL1 
B( 1 , 2)=RIN2/RL2 
B( 1 , 3 )=RIN3/RL3 
RETURN 
END 


£  -j-  y*  y.  y.  y.  y-  y.  y.  y-y  .  Vr  V?  V-  Vr  Vr  Vr  Vr  V.  V?  Vr  Vc  “V  Vr  Vr  Vc 

SUBROUTINE  R2IN(RIN1,RIN2,RIN3) 


COMPLEX  RIN1,RIN2 ,RIN3 

CALL  F2DAT(DB111 ,ANG111 ,DB112 ,ANG112,DB113 ,ANG113) 

PI=4. *ATAN( 1.  ) 

Vlll=10. **( -DB111/20.  ) 

V112=10.  **( -DB112/20.  ) 

V113=10.  **( -DB113/20.  ) 

RAD111=PI/180.  *ANG111 
RAD112=PI/180.*ANG112 
RAD113-PI/ 180. *ANG113 

QRE 1 1 1=V 1 1 1*CQS( RAD  111) 

QRE112=V112*C0S(RAD112) 

QRE113=V113*C0S(RAD113) 

QIM111=V111*SIN(RAD111) 

QIM112=V112*SIN(RAD112) 

QIM113=V113*SIN(RAD113) 

RIN1=CMPLX(QRE111 ,QIM111) 

RIN2=CMPLX( QRE 1 12 , QIM1 12 ) 

RIN3=CMPLX(QRE113 ,QIM113) 

RETURN 

END 


Q  **********-'f******Vf******Vf*****************Hr*^M'f*********yrilr 

SUBROUTINE  F2DAT(DB111 , ANG111 ,DB1 12 , ANG112 ,DB1 13 , ANG113) 

C  ***yr**yryr*****yr*Vf**yr*')’rVfVryf******'>'.-****Vf*****yf**yr*****yr****Vr* 

INPUT  DATA  FILE  IS  DISK  10 
(For  FIXTURE  B) 

READ( 10 , 20)DB 111 ,ANG11 1 , DB112 , ANG112 ,DB113 ,ANG1 13 
20  FORMAT(6(F7. 2) ) 


RETURN 

END 


C  ***************** 

SUBROUTINE  CPX(A) 

C  ***************** 


INTEGER  I PC  3) 

COMPLEX  A( 3 , 3) ,W( 3) ,DET( 2) 

CALL  CG£FA( A , 3 , 3 , IP , INFO) 

CALL  CGEDI(A,3,3, IP,DET,W,01) 

RETURN 

END 


C  ******************************* 
SUBROUTINE  RLOAD(RLl ,RL2 ,RL3 ,F) 

COMPLEX  RL1 , RL2 , RL3 , ZO , ZL1 , ZL2 , ZL3 
REAL  L1,L2,L3 
PI=4.*ATAN(1.  ) 

C=3. *10**8 
Z0S=68. 2 
ERE=2.  833 
ERES=(ERE)**0.  5 
L1=0. 04 
L2=0. 03 
L3=0. 02 

A1=C0S( 2.  *PI*F*L1*ERES/C) 
B1=SIN(2.*PI*F*L1*ERES/C) 

A2=C0S( 2.  *pi*F*L2*ERES/C) 

B2=SIN( 2. *PI*F*L2*ERES/C) 

A3=C0S(2.  *PI*F*L3*ERES/C) 

B3=SIN( 2.  *PI*F*L3*ERES/C) 

R1=-Z0S*( Al/Bl) 

R2=-Z0S*( A2/B2) 

R3=-Z0S*( A3/B3) 

ZL1=CMPLX( 0.  ,R1) 

ZL2=CMPLX( 0.  ,R2) 

ZL3=CMPLX(0.  ,R3) 

ZO=CMPLX( 50.  ,0.  ) 
RL1=(ZL1-Z0)/(ZL1+Z0) 

RL2=( ZL2 - ZO ) / ( ZL2+Z0 ) 

RL3=( ZL3-Z0) / ( ZL3+Z0) 

RETURN 

END 


******************************* 
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THE  INPUT  DATA  FILE  FOR  DISK  8 
3  DIFFERENT  LOAD  FOR  FIXTURE  A 

-V  **  *  *  *  *  ***  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 


3.  07 

22.  7 

2.  28 

98.  2 

2.  14 

175.  1 

3.  82 

-63.  1 

2.  89 

29.  5 

2.  21 

114.  5 

3.  79 

-153. 9 

3.  77 

-49.  2 

2.  57 

55.  4 

3.  46 

126.  4 

4.  12 

-138. 0 

3.  25 

-5.  1 

3.  20 

60.  5 

3.  75 

136.  6 

3.  86 

-74.  3 

3.  33 

-1.  7 

3.  31 

65.  1 

4.  44 

-156. 0 

4.  04 

-69.  9 

3.  40 

4.  9 

4.  42 

119.  8 

5.  39 

-152.  8 

3.  60 

-54.  2 

3.  98 

46.  6 

6.  17 

107.  6 

4.  55 

-119. 0 

3.  89 

-14.  2 

5.  32 

17.  1 

5.  78 

159.  8 

3.  89 

-66.  5 

4.  65 

-51.  0 

6.  39 

61.  8 

4.  15 

-120. 7 

3.  97 

-111.  9 

5.  76 

-25.  5 

5.02 

175.  9 

4.  89 

-175. 9 

5.  70 

-89.  0 

8.  22 

99.  0 

7.  73 

104.  5 

4.  69 

-134. 2 

8.  95 

9.  5 

7.  36 

-17.  1 

3.  91 

173.  7 

6.  37 

-79.  5 

6.  70 

-121.  4 

4.  48 

117.  4 

4.  35 

-147. 2 

4.  59 

170.  4 

6.  84 

47.  9 

3.  49 

157.  8 

4.  12 

114.  2 

8.  48 

-61.  4 

3.  37 

108.  4 

4.  88 

58.  9 

7.  57 

-174. 8 

3.  90 

58.  8 

7.  91 

-9.  1 

4.  92 

115.  1 

5.  55 

1.  4 

8.  42 

-137. 9 

3.  99 

62.  0 

6.  59 

-82.  6 

******************************* 
THE  INPUT  DATA  FILE  FOR  DISK  9 
BOTH  FIXTURE  AND  51  OHM  RES. 

ft************************'/?'***** 


11.  14 

-9.  3 

4.  60 

51.  0 

4.  40 

51.  4 

6.  45 

-150. 2 

11.  38 

-81.  2 

4.  71 

-1.  5 

4.  64 

0.  2 

6.  27 

173.  1 

10.  18 

-154. 3 

4.  80 

-53.  9 

4.  90 

-52.  2 

6.  05 

137.  6 

8.  47 

143.  1 

5.  35 

-106.  1 

5.  36 

-103. 3 

5.  86 

103.  2 

7.  13 

89.  7 

5.  84 

-152. 2 

5.  87 

-150.5 

5.  51 

70.  3 

6.  35 

41.  0 

6.  25 

159.  0 

6.  30 

160.  2 

4.  99 

37.  0 

6.  21 

-5.  3 

6.  99 

112.  1 

6.  99 

113.  1 

4.  68 

3.  0 

6.  51 

-51.  8 

7.  49 

66.  6 

7.  54 

68.  4 

4.  48 

-30.  7 

7.  33 

-102. 0 

7.  25 

19.4 

7.  32 

22.  0 

4.  43 

-65.  6 

8.  92 

-163. 3 

6.  65 

-26.  9 

7.  04 

-25.  9 

4.  94 

-101.  3 

11.  29 

112.  5 

5.  81 

-79.  3 

7.  10 

-73.  4 

5.  75 

-133.  5 

12.  01 

12.  6 

6.  31 

-133.  3 

6.  46 

-121.  3 

6.  18 

-165. 5 

11.  59 

-61.  5 

7.  25 

179.  3 

7.  30 

174.  6 

6.  91 

160.  7 

9.  15 

-93.  7 

6.  90 

135.  8 

6.  45 

133.  3 

7.  83 

125.  8 

6.  27 

-142.5 

6.  56 

85.  4 

6.  35 

83.  2 

9.  07 

87.  1 

5.  53 

169.  5 

6.  25 

34.  5 

6.  00 

34.  1 

1.  78 

41.  0 

6.  07 

126.  3 

5.  85 

-17.  5 

5.  61 

-14.  6 

12.  72 

-22.  9 

7.  54 

87.  7 

5.  50 

-69.  5 

5.  30 

-64.  8 

12.  53 

-104. 0 

9.91 

51.8 

5.  39 

-121.  9 

5.  34 

-118.  0 

10.  00 

-174. 9 

10.  46 

23.  7 

5.  66 

-176. 7 

5.  70 

-171.  1 

7.  82 

131.  2 

11.  95 

44.  3 

6.  05 

131.  4 

6.  10 

134.  9 

6.  53 

87.  5 

w 


vJ 


•»* 


•y 


m 


5 

iS'l; 


i 


$ 


m 


u . 

w: 


:5 


*;vi 


W<l 


89 


“•ft 


j 


SI 


I 


:a 


I 


?S 

ri 


I 


i 


I 


V»f 

>'»*3 


*iV******->V*Vf****1V**iWr>V******TWc** 

THE  INPUT  DATA  FILE  FOR  DISK  10 
3  DIFFERENT  LOAD  FOR  FIXTURE  B 

Vr**********^V**-***Vf***AVrVf****** 


1.  97 

-117. 8 

2.  32 

-26.  9 

1.41 

78.  5 

1.  89 

-176. 2 

2.  46 

-91.  9 

1.  71 

28.  1 

2.  00 

130.  2 

2.  44 

-155. 5 

2.  16 

-26.  2 

2.  27 

81.  9 

2.  51 

147.  3 

2.  48 

-84.  3 

2.  40 

34.  8 

2.41 

96.  6 

2.  69 

-145. 9 

2.  62 

-13.  9 

2.  25 

50.  1 

2.  85 

152.  4 

3.  04 

-67.  3 

2.  35 

6.  9 

2.  85 

94.  8 

3.  87 

-130. 3 

2.  28 

-35.  9 

2.  75 

43.  2 

5.  05 

151.  6 

2.  51 

-78.  5 

2.  77 

-2.  9 

5.  41 

64.  0 

2.  69 

-123.  9 

2.  68 

-44.  3 

4.  91 

-11.  7 

3.  25 

-175.  1 

2.  68 

-82.  4 

3.  87 

-67.  6 

4.  27 

126.  2 

1.  93 

-121.  1 

2.  90 

-118. 9 

5.  24 

56.  1 

1.  91 

-162.  1 

2.  55 

-166. 3 

5.  38 

-20.  6 

2.  02 

157.  0 

2.  68 

147.  0 

4.  49 

-90.  2 

2.42 

113.  2 

3.  33 

96.  3 

3.  49 

-147. 9 

3.  21 

64.  6 

4.  95 

36.  4 

2.  76 

163.  0 

4.  49 

7.5 

7.  33 

-43.  4 

2.  43 

118.  9 

5.  70 

-60.  8 

7.  51 

-139.  7 

2.  36 

77.5 

5.  80 

-134. 9 

5.  35 

143.  6 

2.  55 

36.  2 

4.  89 

158.  8 

3.  67 

87.  9 

3.  32 

-8.  6 

3.  72 

105.  1 

APPENDIX  F.  FORTRAN  PROGRAM  FOR  DE-EMBEDDING  THE 

FILTER 


DIMENSION  X( 21) ,MDB11( 21) , PSE 11(21) 

MDB12(21) ,PSE 12(21) 

MDB21(21) ,PSE21(21) 

MDB 22(21) ,PSE22(21) 

REAL  MDB11 ,MDB12 ,MDB21 ,MDB22 
MAG11 ,MAG12 ,MAG21 ,MAG22 
IM11 , IM12 , IM21 , IM22 
COMPLEX  Sll  ,S12 ,S21 , S22 
All , A12 , A21 ,A22 
B11,B12,B21,B22,W 
Dll ,D12 ,D21 ,D22 
Cll  ,C12,C21,C22 
Fll ,F12 ,F21 ,F22 
Mil ,M12 ,M21 ,M22 
H11,H12,H21,H22 

PI=4. *ATAN( 1.  ) 

DO  10  1=1,21 

INPUT  DATA  FILE  IS  DISK8 
(FIXTURE  1  ONLY) 

READ( 8 ,20)DB11 ,DB12,DB21 ,DB22, ANG11 , ANG12.ANG21 ,ANG22 
FORMAT( 8(F7.  2) ) 

VI 1=10.  **( -DB 11/20.  ) 

V 12=10.  **( -DB12/20.  ) 

V21=10.  **( -DB21/20.  ) 

V22=10. **( -DB22/20.  ) 

RAD11=PI/180.  *ANG11 
RAD12=PI/180.*ANG12 
RAD21=PI/ 180.  *ANG21 
RAD22=PI/180.*ANG22 

QRE 1 1=V 1 l*COS ( RADII) 

QRE 1 2=V 1 2*C0S ( RAD 12) 

QRE21=V21*C0S(RAD21) 

QRE22=V22*COS(RAD22) 

QIM11=V11*SIN( RADII) 

QIM12=V12*SIN(RAD12) 

QIM21=V21*SIN(RAD21) 

QIM22=V22*SIN(RAD22) 

A1 1=CMPLX( QRE 1 1 , QIM1 1 ) 

A12=CMPLX( QRE 12 , QIM12 ) 


A2 1=CMPLX( QRE2 1 , QIM2 1 ) 

A22=CMPLX(QRE22 ,QIM22) 

C  "S"PARAMETERS  TO  "T"PARAMETERS 

B11=-(A11*A22-A12*A21)/A21 
B12=  A11/A21 
B21=-A22/A21 
B22=l/A21 

C  INVERSE  OF  FIXTURE  1  "T" PARAMETER 

W=  1/(B11*B22-B21*B12) 

C11=W*B22 
C12=-W*B21 
C21=  -W*B12 
C22=  W*B11 

CALL  OF  ’ T' PARAMETERS  OF  BOTH  FIXTURE  +  FILTER 

CALL  TMFS(D11 ,D12 ,D21 ,D22 , I) 

C  (INVERSE  OF  FIXTURE  1)  *  (BOTH  FIXTURE+DEVICE) 

H11=C11*D11+C21*D12 
H12=C12*D11+C22*D12 
H2 1=C 1 1*D  2 1+C  2 1*D  2  2 
H22=C12*D21+C22*D22 

CALL  OF  INVERSE  'T' PARAMETER  OF  FIXTURE  2 

CALL  FIX2(F11,F12,F21,F22) 

C  CASCADE  OF  FIXTURE  1,  FILTER  AND  FIXTURE  2 

M 1 1=H 1 1*F 1 1+H2 1*F 1 2 
M 1 2=H 1 2*F 1 1+H2  2*F 1 2 
M2 1=H1 1*F2 1+H2 1*F22 
M22=H12*F21+H22*F22 

C  "T" PARAMETERS  TO  "S"PARAMETERS 

S11=M12/M22 

S12*  M11-(M12*M21)/M22 
S2I=  1/M22 
S22=  -M21/M22 

QIM1=AIMAG(S11) 

QRE1=REAL(S11) 

QIM2=AIMAG(S12) 

QRE2=REAL(S12) 

QIM3=AIMAG( S21) 

QRE3=REAL(S21) 

QIM4=AIMAG( S22) 

QRE4=REAL(S22) 


S2 


DE -EMBEDDED  OUTPUT 


MAGI 1=CABS( S 1 1 ) 

MAG12=CABS(  S12) 

MAG21=CABS(S21) 

MAG22=CABS(  S22) 

MDB 1 1( I )=ABS( 2O*ALOG10( MAG11 ) ) 

MDB12(  I  )=ABS(  20,'rALOG10( MAG  12) ) 

MDE21( I)=ABS( 20*ALOG10(MAG21) ) 

MD322C I)=ABS( 20*AL0G10(MAG22) ) 

PSE11( I )=ATAN2(QIM1 ,QRE1)*180.  /PI 
PSE12( I )=ATAN2( QIM2 ,QRE2)*180.  /PI 
PSE2I( I)=ATAN2(QIM3 ,QRE3)*160.  /PI 
PSE22( I)=ATAN2(QIM4,QRE4)*180.  /PI 

N=5680 

X( I)=N+320.  *FLOAT( I) 

WRITE( 6, 100)  X( I) ,MDB11( I) ,MDB12( I) ,MDB21( I) ,MDB22( I) , 
1  PSE11( I) ,PSE12( I) ,PSE21( I) ,PSE22( I) 

FORMAT! '  '  ,2X,F6.  0,8(1X,F7. 2)) 

CONTINUE 

END 


***iWt*****iV****iWt  *********************** 


SUBROUTINE  TMFS  !D11,D12,D21,D22,I) 
*******  ******************************** 


COMPLEX  A11,A12,A21,A22 
COMPLEX  Dll ,D12 ,D21 ,D22 
PI=4.*ATAN!1.  ) 

INPUT  DATA  FILE  IS  DISK  9 
! BOTH  FIXTURE  +  DEVICE  ) 

READ! 9 , 21)DB11 ,DB12 ,DB21 ,DB22 ,ANG11 ,ANG12 ,ANG21 , ANG22 
FORMAT! 8!F7. 2)) 

Vll=10. **! -DB11/20. ) 

V12=10. **( -DB12/20.  ) 

V21=10.  **( -DB21/20.  ) 

V22=10.  **( -DB22/20.  ) 

RAD11=PI/180.*ANG11 
RAD12=PI/180.*ANG12 
RAD21=PI/180.*ANG21 
RAD22=PI/180.  *ANG22 

QRE11=V11*C0S! RADII) 

QRE 1 2=V 1 2*COS ! RAD 1 2 ) 

QRE2 1=V2 l*COS! RAD2 1 ) 

QRE22=V22*COS!RAD22) 

QIM11=V11*SIN( RADII) 

QIM12=V12*SIN!RAD12) 

QIM2 1=V2 1*SIN( RAD2 1 ) 

QIM22=V22*SIN( RAD22) 


A11=CMPLX(QRE11 ,QIM11) 
A12=CMPLX(QRE12,QIM12) 

A2 1=CMPLX( QRE21 , QIM2 1 ) 
A22=CMPLX( QRE22 ,QIM22) 

D11=-(A11*A22-A12*A21)/A21 

D12=  A11/A21 

D21=-A22/A21 

D22-1/A21 

RETURN 

END 


******************************** 
SUBROUTINE  FIX2(F11 ,F12 ,F21 ,F22) 

'5V-5WoV->V*Vr*->V**,:V**')Wc,*^nWr)V')'r,)Wr^nVyrVpVr^nSnWf 


COMPLEX  All ,A12 ,A21 ,A22 
COMPLEX  B11,B12,B21,B22,W 
COMPLEX  C11,C12,C21,C22 
COMPLEX  Fll ,F12 ,F21 ,F22 
PI=4.  *ATAN( 1.  ) 

INPUT  DATA  FILE  IS  DISK  10 
(FIXTURE  2  ONLY) 

READ( 10 , 20)DB11 ,DB12 ,DB21 ,DB22,ANG11 ,ANG12 , ANG21 , ANG22 
FORMAT( 8(F7. 2)) 

Vll=10. **( -DB11/20. ) 

V12=10. **( -DB12/20.  ) 

V21=10.  **( -DB21/20.  ) 

V22=10.  **( -DB22/20.  ) 

RAD 11=P 1/180.  *ANG11 
RAD12*PI/180.*ANG12 
RAD21=PI/180.*ANG21 
RAD22=PI/ 180.  *ANG22 

QRE11=V11*C0S( RADII) 

QRE 1 2=V 1 2*COS( RAD 12 ) 

QRE2 1=V2 l*COS( RAD2 1 ) 

QRE22=V22*COS(RAD22) 

QIM11=V11*SIN( RADII) 

QIM12=V12*SIN(RAD12) 

QIM21=V21*SIN(RAD21) 

QIM22=V22*SIN(RAD22) 

A11=CMPLX(QRE11,QIM11) 

A12=CMPLX(QRE12 ,QIM12) 

A21=CMPLX( QRE21 ,QIM21) 

A22=CMPLX(QRE22,QIM22) 

B 1 1=- ( A 1 1*A22-A12*A2 1 ) / A2 1 
B12=  A11/A21 


B21=-A22/A21 

B22=l/A21 


V=  1/ ( B 1 1*B22-B21*B 12) 

F11=W*B22 

F12=-W*B21 

F2 1=  -W*B12 

F22=  W*B 1 1 

RETURN 

END 


The  Input  Data  File  for  Disk8 

•jWoWrVc^VVoVVoWfVfVnV^Wf^VVfVrVcVfVcVc^V^WnWcVriV 


14.  48 

.  19 

0.  10 

16.  22 

116.  0 

34.  2 

34.  5 

125.  9 

19.  83 

.  20 

0.  11 

14.  01 

154.  9 

-38.  8 

-39.  5 

148.  9 

22.  21 

.  15 

0.  11 

18.  46 

146.  2 

-113.  3 

-114. 3 

105.  1 

24.  81 

.  10 

0.  14 

19.  94 

142.  8 

170.  8 

170.  2 

139.  2 

14.  41 

.  20 

0.  15 

12.  81 

171.  7 

94.  6 

95.  2 

140.  2 

14.  49 

.  09 

0.  11 

14.  74 

133.  1 

21.  3 

21.  0 

104.  3 

22.  80 

.  05 

0.  09 

19.  44 

80.  8 

-54.  8 

-54.  3 

82.  3 

23.  32 

.  10 

0.  10 

19.  30 

179.  7 

-129.  7 

-129. 3 

78.  0 

17.  71 

.  10 

0.  15 

30.  65 

160.  7 

154.  9 

155.  7 

-16.  9 

18.  60 

.  15 

0.  11 

18.  74 

144.  4 

79.  2 

79.  8 

-162. 4 

14.  23 

.  30 

0.  30 

13.  76 

88.  8 

3.  7 

4.5 

99.  6 

15.  56 

.  26 

0.  30 

16.  14 

54.  1 

-70.  3 

-69.4 

43.  3 

19.  00 

.  17 

0.  27 

20.  31 

-19.  0 

-145. 9 

-144. 3 

-102. 8 

12.  90 

.  10 

0.  20 

27.  01 

-58.  7 

138.  1 

140.  7 

-133. 6 

17.  20 

.25 

0.  20 

28.  54 

-103. 0 

64.  4 

65.  6 

25.  3 

17.  25 

.  19 

0.  24 

20.  37 

-97.  6 

-10.  6 

-9.  5 

-103. 5 

13.  41 

.  20 

0.  25 

15.  86 

-147. 4 

-86.  7 

-85.  2 

170.  9 

16.  90 

.  31 

0.  33 

21.  53 

-169.6 

-160. 3 

-158. 4 

164.  2 

16.  80 

.  25 

0.  20 

25.  37 

-96.  5 

125.  6 

126.  5 

-133.  1 

25.  28 

.  22 

0.  23 

23.  05 

-134. 6 

50.  0 

50.  4 

51.  9 

20.  14 

.  25 

0.  25 

26.  70 

-102. 3 

-27.  9 

-25.  8 

-117. 8 

S5 


iV  -jV  Vr  -k  ir  Ve  ->V  V<  Vr  Vr  iV  tV  Vc  f r  Vr  Vc  t  V  ■>  V  ■>’{■  Vr  -)V  7 Wr  *  Vr  7V  7 V  7V  ->V  Vr 


The  Input  Data  File  for  Disk9 

^VVfVnV^WrsWf’jV^rjV'jWrjViWoV^WoVVfVrjVtWftWfjWoV 


.  84 

87.  81 

88.  98 

.  72 

-150. 6 

104.  6 

134.  3 

54.  2 

.  86 

88.  94 

90.  45 

.  91 

-9.  1  - 

142.  0 

-137. 8 

-113.  7 

.  96 

78.  95 

77.  70 

.  90 

136.  6 

23.  7 

26.  9 

91.  0 

.  75 

71.  67 

71.  43 

1.  03 

-83. 1  - 

163.  8 

-164. 0 

-86.  6 

.  98 

52.  79 

52.  97 

.  62 

45.  7 

-12.  8 

-10.  2 

105.  1 

.  75 

31.  40 

32.  85 

.  83 

150.  7 

128.  9 

116.  1 
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